This Page Is Inserted by IFW Operations 
and is not a part of the Official Record 

BEST AVAILABLE IMAGES 

Defective images within this document are accurate representations of 
the original documents submitted by the applicant. 

Defects in the images may include (but are not limited to): 

• BLACK BORDERS 

• TEXT CUT OFF AT TOP, BOTTOM OR SIDES 

• FADED TEXT 

• ILLEGIBLE TEXT 

• SKEWED/SLANTED IMAGES 

• COLORED PHOTOS 

• BLACK OR VERY BLACK AND WHITE DARK PHOTOS 

• GRAY SCALE DOCUMENTS 

IMAGES ARE BEST AVAILABLE COPY. 



As rescanning documents will not correct images, 
please do not report the images to the 
Image Problem Mailbox. 



REMARKS/ARGUMENTS 



Claims 1-27 are pending. Claims 13, 14 and 27 have been withdrawn. 

The amended first sentence of the specification properly claims priority to the relevant 
patent applications, thereby obviating the objection by the Examiner. 

Applicants gratefully thank the Examiner for her courtesy in conducting an 
unscheduled telephonic interview on February 27, 2004, at which applicants' representative 
and the Examiner discussed, among other issues, that the instant claims are drawn to a 
screening method for determining a genotype associated with bipolar affective disorder, not 
to a genotype which is identified by such a screening method. 

Rejections under 35 USC §112 

The inventors of this application have identified six regions within the human genome 
that are associated with resistance to (protection from), or susceptibility to, familiar bipolar 
affective disorder (BPAD). In each of these regions, the inventors have identified at least 
one, and generally considerably more, markers that exhibit a highly significant linkage to 
resistance to, or susceptibility to, BPAD. These regions represent a relatively limited region 
(generally between about 10 and about 40 cM, as is discussed in more detail below) in which 
one can, e.g. , screen members of a family for these or for further markers that are correlated 
with BPAD. That is, they provide a region in which one can determine one or more markers 
(a genotype) that are linked, in a particular family of interest, to susceptibility to, or protection 
from, BPAD. Alternatively, the regions identified in the present application provide the basis 
for narrowing the size of the region that is linked to BPAD, or for strengthening the statistical 
correlations shown in this application. The above types of analyses can be accomplished by 
screening additional markers; by using the disclosed markers, e.g., the markers recited in 
claims 5, 7, 8, 10 and/or 25; and/or by studying larger and/or different cohorts of subjects. 
The methods of the invention allow the formulation of sufficient detail such that 
determination of an allele(s) associated with increased resistance to BPAD may be 
determined (see, e.g., page 5, lines 25-28 of the specification). 

Thus, the instant claims are directed to screening methods , e.g., a method for 
determining a genotype associated with increased resistance to (protection from) familial 
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bipolar affective disorder (BP AD) ... (independent claim 1); a method for determining the 
contribution of a chromosomal region to the presence of resistance to BPAD. . . (independent 
claim 23); or a method for determining a genotype associated with increased susceptibility to 
familial bipolar affective disorder (BPAD) . . . (independent claim 24). In the method of claim 
1 , for example, a genotype of a (at least one) family member is determined; the BPAD status 
of the family member is determined; the genotype is compared with the BPAD status; and 
from this comparison, it is determined if the genotype is associated with increased resistance 
to BPAD. 

The enablement and the written description rejections in the Office Action appear to 
be directed to claims that recite a genotype which is associated with increased or decreased 
resistance to bipolar disorder (e.g., the result of a screening method for determining a 
genotype associated with increased or decreased resistance to BPAD). This is not what the 
instant claims recite. As discussed in the Reply filed on July 22, 2003, and again below, the 
specification provides both enablement and written description for performing the claimed 
screening methods. 

With regard to the identification of the six regions, the Reply of July 22, 2003 
reviewed in detail that at least one marker within each of the regions is significantly 
correlated with protection from or susceptibility to BPAD. For example, the following 
markers exhibit statistically significant linkage to susceptibility to BPAD: D6S7 on 
chromosome 6; D13S1 on chromosome 13; and D15S45 on chromosome 15; and the 
following markers exhibit statistically significant linkage to resistance to (protection from) 
BPAD: D4S2949 on chromosome arm 4p; D4S397 on chromosome arm 4q; and D11S133 
and Dll S29 on chromosome 11. That is, D6S7, D13S1 and D15S45 are each linked to 
susceptibility to BPAD with a SIBPAL p value of 0.0003 (see, e.g., Figures 3, 4 and 5, 
respectively); D4S2949 is linked to resistance to (protection from) BPAD with a SIBPAL p 
value of < 5xl0" 5 (see, e.g., Figure 6 and page 44, lines 11-13); D4S2949 is linked to 
resistance to BPAD with a SIBPAL p value of O.0001 (see, e.g, Figure 7); and Dll SI 33 
and Dll S29 are linked to resistance to BPAD with a nominal p value of < 5xl0" 5 (see, e.g., 
page 44, lines 15-16). One of skill in the art would know that a SIBAL p value (statistical 
significance of the genetic linkage between markers based on sib-pair analysis) of < 0.001 is 
statistically significant under these circumstances. See also the specification at page 34, lines 
13-15 ("those markers having a test statistic of p O.001 in any one analysis type . . . showed 
evidence for linkage"). The Examiner has presented no reasons or evidence to doubt the 
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assertion that these markers are highly correlated with BP AD. 

The demonstration of these six strongly correlated markers would, in itself, be 
sufficient to identify six regions for further genotypic screening; each region extends for 
about lOcM on either side of one of those markers. It is well-known to those of skill in the 
art that one would reasonably expect that markers located at about lOcM (or, in some cases, 
more) from a marker of interest would be linked to that marker. For an example of current 
practice in the field, see the attached reference: Bosse et al (2003) Circulation 107, 2361- 
2368. The Examinees attention is drawn in particular to Figure 1, which shows linkage data 
for a number of markers. On Chr 2, the linked markers extend from about cM position 225 ± 
20cM, for a span of about 40cM; on Chr 4, the linked markers extend from about cM position 
45 ± 20cM, for a span of about 40cM; on Chr 14, the linked markers extend from about cM 
position 80 ± lOcM, for a span of about 20cM; and on Chr 17, the linked markers extend 
from about cM position 60 ± lOcM, for a span of about 20cM. Clearly, markers that extend 
for about lOcM (or more) on either side of a marker of interest are linked to that marker. 

The present application goes beyond this initial identification of six regions of 
interest, by demonstrating numerous additional markers within each of the regions that 
exhibit strong correlation to BPAD. These data are shown both in the Figures and in the 
Tables in the specification. In general, the Figures present selected data points from the larger 
number of data points presented in the Tables. A summary of data identifying some of these 
additional markers is presented in Appendix A. In response to the request by the Examiner 
on page 4 of the Office Action, the Tables in which these data appear are explained in detail 
in the Appendix. Briefly, on chromosome 6, at least 4 markers have been identified within 
the ~38cM region between 6S344 and D6S589 that are correlated with bipolar disease in a 
statistically significant manner. On chromosome 13, at least 4 markers have been identified 
within the claimed ~12cM region between D13S218 and D13S1 71 that are correlated with 
bipolar disease in a statistically significant manner. On chromosome 15, at least 5 markers 
have been identified within the claimed 20cM region between D15S131 and D15S117 that are 
correlated with bipolar disease in a statistically significant manner. On chromosome arm 
4p, at least 19 markers have been identified within the 33.3cM region between D4S431 and 
D4S404 that are correlated with bipolar disease in a statistically significant manner. On 
chromosome arm 4q, at least 7 markers have been identified within the claimed 42.8cM 
region between D4S402 and D4S424 that are correlated with bipolar disease in a statistically 
significant manner. And on chromosome 11, at least 6 markers have been identified within 
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the claimed region between D11S133 and D11S29 that are correlated with bipolar disease in a 
statistically significant manner. 

In the methods recited in the instant claims, members of a family of interest are tested 
for the presence of markers that lie within the indicated regions. As noted above, it would be 
evident to a skilled worker that markers within about 10 cM on either side of the highly 
correlated markers identified in the application are likely to be linked to those markers, 
thereby defining six linkage regions of interest. In fact the regions of interest are even more 
well characterized. For example, for chromosome 6, markers with statistical significance are 
spread across about 30cM; thus, the region of interest extends for at least this length. For 
chromosome 13, markers with statistical significance are spread across about 12cM; thus, the 
region of interest extends for at least this length. The length of the regions of interest for 
chromosomes or chromosome arms 15, 4p, and 4q extend for at least about 20cM, 33cM, or 
45cM, respectively. 

The identification of the six regions of interest is supported throughout the 
specification. The specification is replete with sophisticated analyses of the data, employing 
a variety of statistical techniques, which support the significance of the linkages. The fact 
that most of the data presented in the application have been published in prestigious peer 
review journals confirms the validity of these conclusions. Two such papers are attached for 
the convenience of the Examiner: Ginns et al. (1996) Nature Genetics 12, 431-435; and 
Girmsetal. (1998) Proc Natl Acad Sci USA 95, 15531-15536. The allegation in the Office 
Action that other studies of linkage in BPAD may have been faulty (e.g., as reported in the 
Berrettini reference cited by the Examiner, and as acknowledged in the specification, for 
example at page 2, line 27 to page 3, line 2) does not cast doubt on the accuracy or 
significance of the present findings. 

The fact that several markers within a region of interest may not appear to be 
statistically significant by individual statistical methods does not cast doubt on the conclusion 
that the region is correlated with BAPD. For markers that lie within regions identified as 
showing significant linkage, but which by themselves do not exhibit such linkage, this result 
likely reflects the fact that these latter markers are less informative. This can result, e.g., from 
less diverse alleles, i.e. fewer alleles. Furthermore, different statistical tests (modelings) 
produce different results. For instance, in linkage analyses, one can model recessive and/or 
dominant inheritance, while in allele sharing (SIBPAL, ASP, etc) the results are less 
dependent on a priori knowing the exact modes of inheritance. Since the genetics 
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(inheritance) working for each gene location was not known when this study was initiated, 
multiple statistical methods (e.g., linkage, sibpal, ASP) were performed. These methods are 
described, e.g., on pages 32-34 and 43-44 of the specification. As expected, different results 
were obtained with the different linkage analyses. The variability in values obtained by the 
different statistical methods is gene specific and varies from location to location. For 
example in Table 5, for marker D13S171, the Zmax dominant value of 0.000 appears highly 
significant, whereas the SIBPAL p value of 0.4905 might not appear by itself to be 
statistically significant. In a case like the above, when the nature of the inheritance was not 
known in advance, it is perfectly acceptable to identify the most significant values, and to 
conclude from them that the marker does, indeed, show linkage to BPAD. (The significance 
obtained under assumptions of a specific model can also be used to identify the genetic 
inheritance mechanisms.) 

The Examiner appears to be concerned that the Blackwood reference cited in the 
Office Action (Nature Genetics 12, 427-430, 1996) identifies some markers in chromosome 
arm 4p as being associated with susceptibility to BPAD, whereas the present inventors 
identify markers in this region as being protective for BAPD. This is not a problem. It is not 
uncommon to find alleles within a given chromosome region, or even mutations within a 
given gene, that result in opposite phenotypic effects. For example, Geller et al. (2004) Am J 
Hum Genet. 74, 572-81 report that some mutations in the melanocortin-4 receptor gene 
predispose a subject to obesity, whereas other alleles in this gene have a negative association 
with obesity. This article is attached for the convenience of the Examiner. Furthermore, 
certain isoforms of ApoE (e.g., ApoE 2) are protective for Alzheimer's Disease, whereas 
others (e.g., ApoE 4) are associated with susceptibility to that disease. Three articles that 
discuss such results with the ApoE alleles are attached for the convenience of the Examiner: 
Sevush et al. (2000) Am J Geriatr Psychiatry 8, 254-256; Beyer et al. (2002) NeuroReport 
13 , 1403-1405; and Higgins et al. (1997) Pharmacology Biochemistry and Behavior 56, 675- 
685. The present inventors have recognized that markers on chromosome arm 4p correlate 
with protection from BPAD, an association that was not reported by Blackwood. 

As for the allegation in the Office Action on page 7 that "The ability to screen for a 
wellness allele is even more unpredictable because it is very difficult to distinguish between 
the presence of a protective allele and the absence of a susceptibility allele," the specification 
clearly teaches how to distinguish between these two possibilities. For example, the 
specification states at page 39, lines 9-16 that 
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Importantly, because of the long-term, longitudinal nature of the study, even the 
unaffected, mentally healthy individuals (those without any psychiatric illness) in 
these families have been closely followed, many for a period of years past the age of 
risk for BPAD. Consequently, rather than limit this genome-wide search to 
identifying susceptibility loci for the disease phenotype (BPAD), we tested the 
hypothesis that "protective" alleles may contribute to the absence of psychiatric 
illness (i.e., mental health 'wellness') in unaffected family members in the 'high 
risk' pedigrees. 

and at page 51, lines 13-16: 

Accordingly, an important step in our study which demonstrates that there are 
"protective" alleles was to show that there are "mentally healthy" individuals who 
share marker alleles that should increase the risk of developing BPAD, and yet, in 
the presence of 'protective' alleles these individuals do not manifest BPAD. 

The existence of markers associated with wellness is well-established in the art. 
See, e.g., the specification at page 54, line 16 to page 55, line 4. For example, wellness 
alleles have been reported for viral infections and Alzheimer's disease. Meyer et al 
(1 997) AIDS LL F73-8 (attached) reports that a CCR-5 delta 32 mutant protects against 
HIV-1 disease progression. Beyer et al, supra and Higgins et ai, supra disclose 
protective effects of ApoE2 alleles for Alzheimer's Disease. Again, the fact that the 
well respected journal, Proc. Natl. Acad. Sci, USA, published the above-mentioned 
paper by the present inventors, which shows that the markers on chromosome 4p from 
the present application are linked to a protective effect for BPAD, attests to the 
convincing nature of those data. 

The specification is fully enabling for the recited screening methods: it teaches how to 
obtain samples from family members (e.g., at page 19, lines 12-27 and Example I); how to 
perform the genotypic analysis (e.g., in Example II); how to assess BPAD status (e.g., in 
Examples I and IV) and how to analyze statistically whether the genotype determined is 
associated with increased or decreased resistance for BPAD (e.g., in Example EH). The 
specification teaches a variety of markers that can be used in the claimed methods, and 
teaches how additional markers can be generated (e.g., at page 23, line 27 to page 27, line 
10). The Office Action states at page 7 that "While one could conduct additional 
experimentation to determine whether markers exist within the recited regions on 
chromosome 4 and 1 1 and these newly discovered markers are associated with BPAD, the 
outcome of such research cannot be predicted." That is to some extent correct. The methods 
of the invention are screening methods, directed toward determining whether such markers 
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are, in fact, associated with BPAD. 

Not only is the specification enabling for the claimed invention, but it also provides 
written description . Contrary to the allegation of the Examiner on pages 12-14 of the Office 
Action, the specification provides even more than a "representative number" of markers 
within each of the claimed regions of the chromosomes. As discussed above, only one highly 
correlated marker is required to identify a region for further screening. In fact, the present 
application provides many more markers within each of the six regions identified. For 
example, some suitable markers to test for resistance to BPAD, flanking D4S2949 (on 
chromosome arm 4p), are indicated at page 14, line 31 to page 15, line 4; some suitable 
markers flanking D4S397 (on chromosome 4q) are indicated at page 15, lines 8-9; and some 
suitable markers on chromosome 1 1, in the -20 cM span between D11S133 and D11S29, are 
indicated at page 15, line 13. See also further candidate markers for resistance indicated on 
page 18, lines 5-13. Examples of suitable primers for identifying the above markers are 
indicated on page 22, line 1 to page 23, line 10. The specification also teaches some suitable 
markers to test for susceptibility for BPAD, and methods for identifying those markers. 

In view of the preceding amendments and arguments, the application is believed to be 
in condition for allowance, which action is respectfully requested. If the Examiner would like 
to confer about the issues discussed herein, she is welcome to contact the undersigned 
applicants' representative. 

Should any additional fee be deemed due, please charge such fee to our Deposit 
Account No. 22-0261 and advise us accordingly. 



Respectfully submitted, 



Date: March 3, 2004 




Nancy J. Axelrod, Ph.D. 
Registration No. 44,014 



Venable 
P.O. Box 34385 
Washington, D.C. 20004-1601 
Telephone: (202) 344-4000 
Direct: (202)344-8334 
Telefax: (202)344-8300 
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APPENDIX A 

Table 5 (pages 35-36 of the specification; regions from chromosomes 6, 13 and 15), 
Table 7 (pages 46-47; region on chromosome arm 4p) and Table 8 (pages 47-48; region on 
chromosome arm 4q) present results of some of the linkage analyses presented in the above- 
referenced application. 

In Table 5, several markers from each chromosome region are examined in four 
categories of subjects (which are defined in the footnote "Categories"). The fourth category 
contains the largest number of subjects, and therefore could provide more significant data 
than the other three groups. The data are presented for several statistical methods: linkage 
analysis (using either dominant ("dom") or recessive ("rec") models); SIBPAL p and ASP p. 
As discussed above, some markers appear to be significantly linked to BPAD in one, but not 
another type of statistical analysis. One of skill in the art would know that, if a given marker 
appears to be significantly linked to BPAD by one of the types of analysis, this is sufficient 
evidence that the marker is, in fact, linked. 

Tables 7 and 8 present SIBPAL analysis of 4p and 4q markers, in either a single 
pedigree (110) or in a total of four pedigrees (1 10, 210, 310, 410). Nominal and simulated 
values are presented for each marker. 

Some of the significant linkages are summarized below: 

On chromosome 6, at least 4 markers have been identified within the -38 cM region 
between 6S344 and D6S589 that are correlated with bipolar disease in a statistically 
significant manner. The data are shown in Table 5: 

D6S344, zmax/dom = 0.000 for each of the 4 patient groups studied; 

D6S7, zmax/rec = 0.0003 for the first patient group; SIBPAL p= 0.0000, 0.0003 or 
0.0001 for three of the patient groups; 

D6S28, zmax/rec = 0.003, 0.000, 0.000, or 0.000 for the 4 patient groups; 

D6S89, zmax/rec = 0.001 or 0.000 for two of the patient groups. 

On chromosome 13, at least 4 markers have been identified within the claimed -12 
cM region between D13S218 and D13S171 that are correlated with bipolar disease in a 
statistically significant manner. The data are shown in Table 5 and/or in Figure 4: 

D13S218, Table: SIBPAL p=0.0171 or 0.0175 for two of the patient groups; 

DUDSlrn, Fig. 4, SIBPAL p=0.0003; 

D13DSU, Fig. 4, SIBPAL p=0.0003; 

D13S171, Table: Zmax/dom=0.000 for all 4 patient groups. 
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In addition, D13S221, which is outside the region that is currently claimed, is shown in Table 
5 to have values of zmax/rec = 0.012 and Zmax/dom=0.000. 

On chromosome 15, at least 5 markers have been identified within the claimed 
~20cM region between D15S131 and D15S117 that are correlated with bipolar disease in a 
statistically significant manner. The data are shown in Table 5 and/or Fig. 5: 

D15S36, Fig. 5 and Table 5, SIBPAL p=0.01 14 or 0.01 18 for two of the patient 
groups, Table 5, zmax/dom=0.000 for the 4 th patient group; 

D15S45, Fig. 5 and Table 5, SIBPAL p=0.0003 or 0.0018 for two of the patient 

groups; 

D15S148, Fig. 5 and Table 5, SIBPAL p=0.0217, 0.0219, 0.0123 or 0.0360 for the 4 
patient groups; 

135yc, Fig. 5, SIBPAL p=0.102 

D15S117, Table 5, zmax/dom=0.000 for the 4 th patient group. 

On chromosome arm 4p, at least 19 markers have been identified within the 
~33.3cM region between D4S431 and D4S404 that are correlated with bipolar disease in a 
statistically significant manner. These data are shown in Table 7. In a few cases, significant 
values for markers are also shown Fig. 6: 

Marker Fig. 6 (SIBPAL p) Table 7 (SIBPAL p. nominal, pedigrees 110-410) 



D4S431 

D4S2366 

D4S3007 

D4S394 

D4S2983 

D4S2923 

D4S615 

Afmanl80xa9 

D4S2928 

D4S1582 

D4S107 

D4S3009 

D4S2906 

D4S2949 

D4S1582 



0.0002 
0.0015 
0.0035 



0.0000 
0.0005 



0.0110 
0.0002 

0.0014 (pedigree 110) 

0.0007 

<\x\0 A 

0.0003 

^xlO" 4 

<lxl0 -4 

<5xl0" 3 

0.0032 

<5xl0" 5 

0.0001 

0.0004 

<lxl0- 7 
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D4S2942 



<lxlO^ 



D4S3048 



0.0036 



D4S419 



0.0002 



0.0004 



D4S404 



0.0004 



0.0004 



Also, D4S391, which lies outside of the claimed region, is shown in Table 7 to have a 
SBPALp=0.0001 (pedigree 110). 

On chromosome arm 4q, at least 7 markers have been identified within the claimed 
~42.8 cM region between D4S402 and D4S424 that are correlated with bipolar disease in a 
statistically significant manner. The data are shown in Table 8 and/or in Figure 7: 

Marker Fig. 7 CSIBPAL p) Table 8. (SBPAL p. nominal, pedigree 1 10) 

D4S402 O.0001 

D4S409 0.0001 

D4S422 0.0108 

D4S175 <0.0001 0.0471 

D4S1579 0.0015 

D4S397 <0.0001 3xl0" 7 

D4S424 0.0003 0.0481 

On chromosome 11, at least 6 markers have been identified within the claimed region 
between D11S133 and D11S29 that are correlated with bipolar disease in a statistically 
significant manner. The six markers are D11S394, D11S133, D11S147, CD3D, D11S285 and 
D11S29. The specification teaches on page 15, lines 10-13 that these six markers are 
significant, and indicates explicitly that D11S133 and D11S29 are linked to resistance to 
BP AD with a nominal p value of < 5xlO" 5 (see, e.g., page 44, lines 15-16). The Examiner has 
presented no reason to doubt the statement by applicants that all 6 markers exhibit statistically 
significant correlation. 

DC2-DOCSl-528217vl 
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ABSTRACT Bipotar affective disorder (BPAD; manic- 
depressive illness) is characterized by episodes of mania and/or 
hypomnnia interspersed with periods of depression. Compelling 
evidence supports a significant genetic component in the sus- 
ceptibility to develop BPAD. To date, however, linkage studies 
have attempted only to identify chromosomal loci that cause or 
increase the risk of developing BPAD. To determine whether 
there could be protective alleles that prevent or reduce the risk 
of developing BPAD, similar to what is observed in other genetic 
disorders, we used mental health wellness (absence of any, 
psychiatric disorder) as the phenotype in our genome-wide 
linkage scan of several large muttigeneration Old Order Amish 
pedigrees exhibiting an extremely high incidence of BPAD. \Vc 
have found strong evidence for a locus on chromosome 4p at 
D4S2949 (maximum GENE hunter -PLUS non parametric linkage 
score = 4.05, P = 5,22 x lO -4 ; SIBPAL P tMpiri cai value <3 x 10" 5 ) 
and suggestive evidence for a locus on chromosome 4q sxtD4S397 
(maximum c ENEHUNTER-P L US nonparametric linkage score = 
3.29, P = 2.57 x 10" 3 ; SiBPALP erop i rica ! value <1 x 10 _J ) that are 
linked to mental health wellness. These findings are consistent 
with the hypothesis that certain alleles could prevent or modify 
the clinical manifestations of BPAD and perhaps other related 
affective disorders. 

Bipolar affective disorder (BPAD) afflicts approximately 1% of 
the population and is associated with a high risk of suicide (1). 
Twin, family, and adoption studies have provided strong evidence 
for an important genetic component in the susceptibility to 
develop BPAD (2), but unlike for other common medical ill- 
nesses, robust biological markers have not been identified for 
BPAD, and genetic linkage studies have had to rely on categorical 
diagnoses. Genetic heterogeneity, phenocopies, genotyping er- 
rors, and the complexities of performing and interpreting statis- 
tical analyses may contribute to some of the inconsistences 
observed in these linkage studies (3-12). However, because the 
inheritance of BPAD is probably multifactorial, the possible 
involvement of multiple genetic componenis of small effect 
and /or the occurrence of major allelic effects only in epislasis 
must be considered. We previously have suggested (13), that in 

The publication costs of this article were defrayed in p;irt by parie charge 
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addition to susceptibility alleles, there could be alleles that reduce 
the risk of developing BPAD in a manner similar to that reported 
for other complex genetic disorders. If model-based linkage 
analyses are used, a false-negative linkage finding could result 
when individuals inherit disease susceptibility alleles but do not 
manifest the phenotype because of the presence of protective 
alleles. The inclusion of individuals who inherit susceptibility 
alleles but do not manifest disease because of protective alleles, 
or of individuals who inherit protective alleles but nevertheless 
manifest the disease, also will reduce the power of model-free 
(allele-sharing) analyses. Thus, regardless of whether model- 
based or model-free analyses are used, wellness or protective 
alleles could have a significant impact on linkage analyses. 

Ascertainment of psychiatric disorders and health among 
several large multigenerational Old Order Amish pedigrees cov- 
ers a period of more than 20 years. Throughout this longitudinal 
study, procedures for assessing and diagnosing subjects have 
remained constant and have included a thorough evaluation of all 
bipolar I (BPI) probands and their relatives (14-16). Morbid risk 
analyses have demonstrated a high prevalence of affective disor- 
der among firstndegree relatives of bipolar probands in these 
families ( 1 7). Because of the long-term, longitudinal nature of the 
study, the unaffected, mentally healthy individuals in these fam- 
ilies also were followed, many for a period of years past the age 
of risk for BPAD. Consequently, rather than limit our genome- 
wide search to identifying susceptibility loci for the disease 
phenotype (BPAD), we tested the hypothesis that protective 
alleles may contribute to the absence of psychiatric illness (i.e., 
mental health wellness) in unaffected family members in these 
high-risk pedigrees. Because the mode of inheritance of any 
gene(s) modifying the relative risk for affective disorder is un- 
known (2) wc relied exclusively on model-free linkage analyses. 
We now report strong evidence for linkage of DNA markers on 
chromosome 4p to mental health wellness in relatives at high nsk 
for. but who did not develop, major affective disorder in several 
large multigenerational Old Order Amish pedigrees with an 
extremely high incidence of BPAD. 



Abbreviations: BPAD, bipolar affective disorder; NPL. nonparameirtc 
linkage; BPI. bipolar I; GH-PLUS. oesehunter-plus; IBD. identical 
by descent. , 
*To whom reprint requests should be addressed at: Clinical ^eur.'- 
science lirnnch. National Institute rf Menial U-:*UU. » ,u l k,, " l 7. l , 1 ' 
Is i i(ioi flll-.nir.. 40 Coment Drive. MS«'-UIJ>. Iicihe*d.i. Ml) -U^-' 
J l<>5 c-mail: yumsc'v up.nimlt r.ili 
5|'.Si I. ami K A.I 1 , omtribdied cqti.iUy hi this wnrk. 
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MATERIALS AND METHODS 

Diagnostic Assessment. Our genctic-epidcmiologic study of 
BPAD among the Old Order Amish in southeastern Pennsylvania 
has been described in detail (18), including the methods for 
ascertainment and diagnostic documentation with informed con- 
sent (15, 19). Diagnoses were made, by using strict research 
diagnostic criteria (16), by a five-member psychiatric review board 
whose members were blind to pedigree membership, diagnostic 
opinions, treatment data from abstracted medical records, and 
genetic marker status. Twenty-five nuclear families were devel- 
oped based on one or more confirmed cases of B PI, which formed 
the structure of pedigrees 110, 210, and 310. Our present report 
uses all of these earlier subjects plus additional expansions, 
especially in pedigree 410, so that our overall study now contains 
346 samples, including those from 50 BPI individuals. The 
unaffected individuals (mentally well or healthy) arc those for 
whom all SADS-L (schedule for affective disorders and schizo- 
phrenia-lifetime version) interview responses were negative, and 
no contradictory reports were given by family informants. Any 
individuals for whom some symptomatology was identified, even 
if it did not meet criteria for a formal diagnosis, were labeled as 
unknowns in our linkage analyses. 

The method used for this longitudinal study is ethnographic 
and culturally appropriate to the field setting. Several members of 
each nuclear family with a BPI proband (BPI nuclear family) are 
seen annually. Those diagnosed with BPI or other major affective 
disorder undergo a yearly course-of-illncss update. Parents of 
each BPI patient are visited regularly and have proven to be 
accurate informants. At least one unaffected sibling (control 
sample) of the married BPI patients has been interviewed yearly 
since 1990 in connection with a prospective study of children at 
risk for BPAD. Hence, at least three members and occasionally 
all members of each BPI nuclear famity have been evaluated 
yearly. Individuals are interviewed anew with the complete 
SADS-L schedule whenever any abnormal mental or emotional 
symptoms are identified by the follow-up mechanisms. The 
long-term, systematic follow-up of the families in our study has 
demonstrated that onset of illness in the Old Order Amish usually 
is reported by multiple informants. 

Patient-Samples. Blood samples were generally collected after 
each first-degree relative (including parents, siblings, and children 
older than age 15) of the BPI probands had been interviewed with 
the complete SADS-L schedule. Samples were obtained with 
written informed consent and coded to maintain confidentiality. 
The phlebotomist was kept blind to pedigree relationships and 
diagnostic status. Lymphoblastoid cell lines, on an average of 
eight members of each nuclear family, were established at the 
Coriell Institute for Medical Research, Camden, N J. and/or the 
Clinical Neuroscience Branch, Intramural Research Program, 
National Institute of Mental Health, Bethesda, MD. The N1GMS 
Human Genetic Mutant Cell Repository catalogue (20) contains 
updated pedigree and diagnostic information for several of the 
Amish pedigrees used in our study. 

Cenotyping. Genomic DNA was obtained from peripheral 
blood samples and/or immortalized lymphoblastoid cell lines ns 
described (9). The order of typed markers on our mapping panels 
was obtained from the genetic location database (21). 

DNA panels for PCR were in a 96-microtiter plate format, and 
aliquoted by using a BioMck robot (Beckman Instruments). The 
DNA samples were PCR-amplified separately with each of the 
DNA markers and then the PCR products were multiplexed, six 
markers per lane, for electrophoresis using GENESCAN software 
on ABI 373 instruments (Applied Biosystems Division, Pcrkin- 
FJmer). The genotypes subsequently were analyzed with GIINO- 
typuk. DNA from several individuals was represented multiple 
times in the gi-ruitypiug panels as cunt nils to evaluate I In- 
consistency of genotypes across different gels. 

Genetic Analysis Software (3S) was used to identity pioblem 
atic marker data, and a utility written in srss (SPSS. Chicago) 
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generated a list of samples that needed to be rerun because of 
inheritance discrepancies or unreadable signals. We maximize the 
useful information by repeating the genotyping/analysis cycles 
until all possible DNA marker genotypes are obtained. Once 
genotyping for a marker was finished, the data were reanalyzed 
with GAS., and observed allelic mutations and other noninher- 
itances were zeroed out in the data file and noted. Histograms 
indicating the marker allele size bins were generated, fastunk 
(22) was used to reanalyze the data before further statistical 
analyses. 

Linkage Analyses. Model-free linkage analyses were con- 
ducted by using the two-point affected sib pair analysis program 
S.A.G.E. sibpal (23) and the multipoint analysis program gene- 
HUNTER-PLUS (GH-PLUS) (24). Because there were a few sib- 
ships with incomplete marker information, marker allele fre- 
quencies were estimated from the entire Old Order Amish family 
data set by using a maximum likelihood method implemented in 
the program mendel/userm13 (25, 26). sibpal was used to 
identify markers showing an excess of alleles shared identical by 
descent (IBD) among unaffected, mentally healthy sib pairs. 
Under the null hypothesis of no linkage between a trait and 
marker, sib pairs would be expected to share on the average 5U * 
of alleles : IBD, but when a trait and marker axe linked 1BU 
sharing will be increased in both affected and unaffected sibpairs. 
Because sibpal assumes marker allele frequencies appropriate 
for random samples, it underestimates the proportion of alleles 
shared IBD by concordant sib pairs when there is linkage. 
Multipoint analyses by using the model-free linkage program 
GH-PLUS produced nonparametric linkage (NPL) scores along 
points at the chromosomal region of interest. Two scoring func- 
tions arc available in GH-PLUS: IBD sharing can be assessed 
among concordant relative pairs (NPL^) or it may beassessed 
among larger groups of concordant relatives (NPUu). Uur anal- 
yses were conducted by using the NPUn statistics as Kruglyak and 
colleagues (24) have demonstrated that the NPUn statistic results 
in a more powerful test than the NPLp^ statistic. 

RESULTS 

First, we analyzed our genome-wide scan dataset looking for 
evidence of chromosomal loci linked to mental health wellness 
(the absence of any psychiatric illness). In these-analyses, only 
mental health wellness, in individuals who were over 45 years ; ol 
age and had a first-degree BPI sibling in their family (^^« 
110, 210, 310, and 410). was the linkage phenotype of interest 
(concordantly unaffected pairs) by using SIBPAL. Of more than 
980 DNA markers, only five markers representing three chro- 
mosome regions had test statistics that were sufficiently ^outlying 
and that were likely to represent significant linkage results. Of £e 
markers on chromosome 4p f D4S2949, which * 
vicinity of the BPAD susceptibility locus reported jVjf^^ 
et tit. (11). had an empirical SIBPAL P value <3 X 10 5 (nominal 
Pvalue<l X 10- 7 ).The marker on chromosome 4q had 

an empirical sibpal P value - 9 X 10" J (nominal /> value - 3 x 
10- 7 ) On chromosome llq. two DNA markers (DJ1SU3 ana 
DUS291 located over an approximately 20<H region, each had 
a nominal P value <5 x 10' 5 (sibpal; simulations were not 
performed). To supplement standard criter.a for assess.n i uie 
significance of our linkage analysis results, we used g'»P" ™ 
techniques (Fig. 1) and the empirical assessment of P values 
(27-29). If each marker assessed in a pairwise linkage analysis is 
unlinked to the trait, then the P values associated with i those 
markers should be uniformly disiributcd. In addition, im. w 
statistics used to generate these P values (for instance f esis i 
the case of sibpal) should Follow ah appropriate distnb ut ion 
plot (generated by usin- proc CHART, SAS) of the test siai sno 
obtained from each pairwise linkage analysis is shown in o- j 
The pint in the inset depicts a line thai should be Itm.. r • 
markers arc unlinked. As seen in 1. there are tmlly * 
statistic values that lifcclv represent false null hypolhcsc. u ■ . 
evidence for significant Imkagc results. In addition, in inc 
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FiC. I. Plot of test statistics obtained from the pair-wise linkage 
results. (It act) A cumulative plot of P-values whose linearity would reflect 
uniformity in /'-values associated with multiple linkage results whose null 
hypotheses were all true (see text). The outlying test statistics and values 
(denoted by arrows) were associated with markers, D4SW7 {i ~ 6.24), 
D4S2949 (/ = 7.79). D4S292S (t = 5.03). D I IS 133 (r = 6.00). and Dl 1S29 
<r = 6.32). 

Fig. I, the small upturned portion of the Rvalue plot near values 
of 1 -- P - 1 represents departures from uniformity and hence 
most likely reflects false null hypotheses. Because of the effort 
required to investigate the significance of these findings and the 
previous evidence for a BPAD-related locus on chromosome 4 
(11), we chose to examine DKA markers on chromosome 4 first 
for linkage to mental health wellness. 

To evaluate the findings on chromosomes 4p and 4q in more 
detail, we genotyped the subpedigrees and nuclear families 
containing at least one sibling with BP1 (Table 1) by using 
additional DNA markers in these interesting regions. Compared 
with our previous report (9) a larger number of individuals was 
included in these analyses (Table 1). In this report, model-free 
linkage analyses using Sibpal and GH-PLUS (24) were per- 
formed by using mental health wellness as the linkage phenotype 
(Tables 2 and 3). In our analyses, individuals having a psychiatric 
diagnosis other than BPI, as well as those having psychiatric 
symptoms but no diagnosis, were classified in the unknown 
category for affected status. In the Amish Study sample of BPI 
patients (n = 50) the mean and median ages of onset (research 
diagnostic criteria) are 24 and 22 years, respectively. Hence, in all 
analyses we used a conservative age cutoff of 45 years to define 
family members with the unaffected wellness phenotype. We also 

Table 1. Old Order Amish subjects inc luded in linkage analysis 

Analysis categories 



Subjects 



Mentally healthy 



Unknowns 



Pedigrees 110. 210. 310.410 
2:25 years old 
2:35 years otd 
2:45 years old 
2:55 years old 
Pedigree 1 10 only 
2:23 years old 
2:35 years old 
245 years old 
=55 years old 



13S 
109 

74 
52 

AS 
37 
31 
23 



S5 
114 
149 
171 

32 
40 
40 
54 



The category of unknowns includes individuals of unknown phenotype. 
individuals with psychiatric diagnoses other than BPI. and individuals who 
are mentally hcalihv but are younger than the particular age cut-ofr used 
in analyses BPI individuals are mil included in the unknown phenotype 
category. In pedigrees 1 10. 210, 310. and 410, 39 people were diagnosed 
with DIM. ci«ht with BPI I. 21 with recurrent depressive disorder, two with 
unipolar deprive disorder, and 15 with : other psychbiric illness. In 
pcdi«ree I U) alone. IS people were diagnosed with HIM. two with [JP1I. 
10 with major depressive disorder, and five with other psyrhianc illness. 
Note- the iiutividuaUurwil in these link;i 3 c analyses represent only ;. subset 
of the entire Amish bipolar pedigrees because only nuclear fanulie:; and 
subpedigrces containing a sibling with I MM were included. 



examined the influence of different age cutoffs for defining well 
individuals and the contribution of different subpedigrees {fam- 
ilies from pedigrees 110, 210. 310, and 410 versus only families 
from pedigree 110) on the test statistics for linkage. Well indi- 
viduals younger than the specified age cutoff were considered to 
have an unknown affected status in (he analyses. 

On chromosome 4p, (he maximum multipoint NPL values 
(GH-PLUS; including only individuals > age 45 years) were 4.05 
(P « 5.22 X 10-<) for pedigree 1 10 and 4.05 (P = 1.84 x 10" J ) 
for all pedigrees, respectively. The maximum multipoint NPL 
values (GH-PLUS; including individuals > age ; 4: > years) for 
markers on chromosome 4q were 3.29 (P = 2.57 x 10 ») for 
pedigree 110 and 2.82 (P = 4.43 X lp for all ped.grees, 
respectively. The GH-PLUS -log w P value as a function of the 
map position at these locations on chromosome 4 is shown in Fig. 
2. sibpal test statistics for markers on chromosomes 4p and 4q are 
shown in Tables 2 and 3. On chromosome 4 the lowest (nominal) 
P values obtained from the SIBPAL test statistics ^"J^.™^" 
D4S2949(*p;P<l X 10' 7 ) and D4S397 (4q;/>« 3 X 10 The 
nrnimum multipoint NPL values (GH-PLUS; including individ- 
uals > age 45 years) for markers on chromosome llq were 2.43 
for pedigree 1 10 and 2.49 for all pedigrees, respectively. 

To obtain empirical lvalues, we simulated genotype data by 
randomly assigning marker alleles to the Sunders and then 
assigning alleles to Their descendants following Mende han inher- 
itance. Allowing for consanguineous matings, the en ire family 
structure was used in marker assignment, thus taking into account 
all relationships between individuals in the dataset. For each 
simulation, after marker assignment, the pedigrees were tr^mcd 
down to that of the nuclear families used in the : linkage analysis 
S.BPAL then was run on the trimmed dataset, and test statist c s for 
concordant and discordant sib pairs were ob / ain f^^5 n ^ 
value is simply estimated as the proportion of "P'^^ J ? 
the simulated statistic is greater than or equal to the observed 
statistic, i.e., the probability that the observed result o 
more extreme would be obtained by * a ^ lo ~-^ 
were conducted for markers on chromosomes 4p and 4q For jsach 
marker, 100,000 replicates were obtained. ^^Tl^Sf 
on chromosome 4p clearly meets the proposed criteria of sign* 
icance for linkage (30). 

DISCUSSION 
If alleles exist that are associated with mental health , wellness 
then the identification of chromosome regions con ta.mng these 
alleles would be enhanced by studying the genetie.ll y nsk. 
mentally healthy members of large, multigenerat, ^.ElK 
like our Old Order Amish families. However, in try ng to identify 
protective or wellness alleles, one must recogn^ ttar there . are 
phenocopies that need to be considered in our ^J^ges 0 ' 1 , 6 
ihe extremely high risk for developing d.sease. "^^"^ 
are undoubtedly well because they do .not inherit any or aU) o 
the requisite suLep.ibi.ity ajelesfor f^^^t 
greatest liability for onset of BPAU in me uiu « _ 
from early teeri through 24 years of age, the 
the well phenotype for individuals who eventual iy w. (develop 
BPAD would be greatest through th,s age penod '« lh " e 
Order Amish families susceptibility alleles for BPAD ^p ob«W> 
.occur in very high frequency. An important step ,n ^ onstra '^ 
that there are protective alleles will be to show ha iherca e 
mentally health? individuals who share marker alleles tha should 
increase the risk of developing BPAD.andyet. .n the P^" ce ° 
protective alleles, these individuals do not manifest BPA1* ■ ho 
many of the markers, fl. an underestimate of the propo t.on ot 
alleles shared IBD in well sibpnirs. mcreases with a =£ 
,c. a more stringent *™*«^^\ffi^MS 
example, with respect to marker IUS29-l9on ip. 1 " " . , 
and I foe age cu.off ,« in« of 25, ^> W 
This finding suggests thai increasing llu. ^ «> r 
ikiics wmc age-related ucll phenocopies. 
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Tabic 2. Results of sidpal Analysis of 4p markers 

_ Pedigree 1 10 Pedigrees 1 10. 210, 310, 410 

/'-value f-valuc 



Marker 


fr (s.e.) 


Nominal 


Simulated 


n (s.e.) 


Nominal 


Simulated 


D4S4I2 


.4749 (.0621) 


.6555 


np 


.5 116 (.0539) 


.4154 


np 


D4S43I 


.5734 (.0441) 


.0523 


np 


.5921 (.03SS) 


.0110 


np 


D4S2366 


.67Sl(.0452) 


.0002 


.0005 


.6024 (.0356) 


.0027 


.0094 


D4S2935 


.5066 (.0218) 


.3825 


np 


.4998 (.01 9S) 


.5043 


np 


D4S3007 


.6233 (.03S6) 


.0014 


.0023 


.5632 (.0337) 


.0330 


.0496 


D4S394 


.6782 (.05 13) 


.0007 


.1)012 


.5955 (.0421) 


.0135 


.0249 


D4S29S3 


.7219 (.04S4) 


<1 X 10~ 4 


np 


.6090 (.0377) 


.0025 


np 


D4S2923 


.6661 (.0446) 


.0003 


np 


.5902 (.0307) 


.0022 


np 


D4S615 


.7161 (.0393) 


<l x 10" J 


np 


.6223 (.0324) 


-.0002 


np 


AfmalS4za9 


.7396 (.0446) 


<! X 10~ 4 


np 


.6220 (.0370) 


.OOOS 


np 


D4S292S 


.7333 (.0257) 


<5 X 10~ 5 


np 


.6369 (.0272) 


<5 x.10" 5 


np 


D4SI605 


.5453 (.025S) 


.0440 


.0472 


.5795 (.0244) 


.0011 


0.0058 


D4SI5S2 


.67S7 (.0616) 


.0032 


.0112 


.6269 (.0557) 


.0139 


.0510 


D4SW7 


.6557 (.0246) 


<5 X 10"* 


.0029 


.6514 (.0243) 


<5 x 10" 5 


.OOSS 


D4S3009 


.7325 (.0552) 


.0001 


np 


.6237 (.0379) 


.OOOS 


np 


D4S2906 


.6460 (.0396) 


.0004 


n P 


.5S53 (.0327) 


.0055 


np 


D4S2949 


.7077 (.0202) 


<1 x 10" 7 


<3 x 10" 5 


.6SSS (.0243) 


<1 x 10" 7 


<3 x 10~ 5 


AfmO$7i S 5 


.5229 (.036S) 


.26S6 


np 


.5114 (.0246) 


.321S 


np 


D4S2944 


.5647 (.0263) 


.0093 


np 


.5428 (.0255) 


.04S3 


np 


D4S403 


.6032 (.0492) 


.0217 


.0233 


-59S9 (.0443) 


.0232 


.0350 


D4S2942 


.7196 (.030$) 


<1 X 10~ J 


np 


.6627 (.0243) 


<1 x 10" J 


np 


D4S29S4 


.5510 (.0396) 


.1032 


np 


.5493 (.0297) 


.0505 


np 


D4S1602 


.6001 (.0561) 


.0412 


np 


.5703 (.03S.1) 


.0356 


np 


D4SI511 


.6242 (.04S9) 


.0077 


np 


.5779 (.0315) 


.0079 


np 


D4S23I1 


.7429 (.0279) 


<5 x 10" 5 


np 


.6327 (.0336) 


.0001 


np 


D4S304S 


.6628 (.0573) 


.0036 


np 


.599S (j0403) 


.007S 


np 


D4S419 


.59$ I (.0270) 


.0004 


.0010 


-5772 (.0319) 


.0100 


.0201 


D4S404 


.67S5 (.04S9) 


.0004 


.0010 


.64 2S (.0470) 


.0020 


.0072 


D4S39I 


.700S (.04S7) 


.0001 


.0003 


.65S5 (.0470) 


.OOOS 


.0035 



fl is the estimated proportion of alleles shared identical by descent, np. simulations not performed. 



It is conceivable that virtually all cases of affective disorder in 
these families are caused by a common set of susceptibility alleles. 
The wellness or protective loci that we tentatively have identified 
could harbor alleles that prevent the manifestation of a BPAD 
spectrum phenotype, which also could include major depressive 
disorder. In our analyses the strongest evidence for protective 
alleles comes from pedigree 110, suggesting that such alleles may 
be more likely in this branch of the family. However, highly 
significant test statistics and multipoint logarithm of odds scores 
(using GH-PLUS) also are observed when pedigrees 110, 210, 
3 10. and 4 10 are used for analyses (Fig. 2 A and B). The decreased 
sharing in proportion of alleles IBD for discordant pairs provides 
further support for the existence of alleles associated with the 
absence of affective disorder (mental health wellness) in these 
families (data not shown). In addition, epistatic interactions 
between alleles also could prevent or del ay onset of an illness such 
as major depressive disorder from developing into BPAD. In- 
deed, as we increase the age-of-risk cutoff for defining the well 
phenotype from 25 to 45 years in our linkage analyses, the number 
of mentally healthy members decreases as expected, yet the 
evidence for linkage increases (data not shown). 

There is some debate on the analysis of sibling pairs as to 
whether the use of inbred sibling pairs results in an increased 
number of false-positives if allele-sharing-based statistical meth- 
ods arc used (31). However, the arguments that (r) inbred sibling 
pairs are likely to share more genes than noninbred sibling pairs 
(i.e., have a kinship factor greater than 0.5) and («) that greater 
regions of the genome would show significant deviations from the 
expected noninbred sibling sharing value of 0.5 are incorrect 
whoa one is merely considering an analysis of sibling pairs 
involving only the transmission of alleles from p.irer.:* to oil- 
spring. The transmission of .illeles I'uim parent * t<< o:upi i;w will 
follow Mendclian ratios, and tins iho null values lor n. l'or 211*0 



sibling allele sharing in any population will be 0,25. 0.50, and 0.25 r 
whenever only parental and sibling genotype information is used. 
However, if the origin of the parental alleles is taken into 
consideration, then there will be greater information about alleles 
shared by sibling pairs from inbred populations. For example, this 
increased information has the potential to resolve ambiguities in 
the sharing of alleles transmitted from homozygous parents, 
because the two copies of the allele in an inbred homozygous 
parent could be IBD. This information also could help resolve 
alleles shared by siblings identical instate into alleles shared IBD, 
showing that alleles transmitted to two offspring from different 
parents may be copies of the same allele because of the related- 
ness of the parents. If genealogy is tnken into account, then the 
increased ability to resolve ambiguities in allele sharing would 
result in greater power in the analysis of inbred sibling pairs (31). 

Ultimately, if inbreeding exists in a population from which 
sibling pairs have been gathered, but one ignores genealogical 
information by merely studying the transmission of alleles from 
parents to offspring, then no increase in false-positive linkage 
results will occur, because Mendel's law applies to inbred as well 
as outbrcd parent-offspring allele transmission studies. On the 
contrary, a decrease in power may result from inbred sibling pair 
analyses because spouses may manifest greater homozygosity and 
therefore provide less informative genotypes for pnrent- 
offspring-based linkage studies. . . 

Genetic mapping of complex disorders with multifoctona 
inheritance could be especially difficult if. in addition to suscep- 
tibility alleles, individuals inherit protective alleles Hint prevent or 
reduce the risk of manifesting the disease phenotype. fcvc 
thou-h model-based linkage analyses that do not allow lor < 
nmltiUMoiial component are of only limited usefulness in n«->- 
araun,i:>no:v I hey :n'e still mod frequently. In these ,,lM [ Ul ^j 1 
t.iUe-ncvaliVe linkage findim: (lypc 2 erior) could refill «■ 
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Marker 




r cuigrcc mi 






Pedigrees 110. 210. 310. 410 




(1 (s.c.) 




F-vatuc 


i\ (s.c.) 


P-vaUtc 




Nominal 


Simulated 


Nominal 


Simulated 


t/TJ _„'l/J 


4fi If, f AJ">1> 


SMS 


np 


4670 f 0305^ 


.35S5 


nn 






00 ">7 


np 


5403 i* 01391 


.0025 


no 
r 




c j i; / f\t l c\ 
j (.v-1 IJJ 


. 1 J 


.1373 


5445 f 0304 1 


.0743 


•0S46 




.J J J j (.UJ- -J 


.1570 


np 


5"5 f 03S1) 


.2780 






.XJ.O (.U.OVJ 


.«IO IV 


np 




.6370 


no 


I^Vj / /J 


.jvuQ (.IJJOo) 






5995 ( 04541 


.0231 


.034S 


D4h422 


.6I9S (.OjOU) 


A ■ AO 


np 


5A<15 f 0VR61 


.0403 


nn 


D4SI576 


.5290 (.0:>09) 


.286 1 


np 


5177 / 01671 


. 1 545 

" .6358 


no 


D4S22Q4 


.4960 (.0446) 


.XO 1 


np 


d<!/;7 / 03SH 

.HOO / ^.uJu 1 ) 


nD 


n if if 7n 


.6206 (.Ojol) 


nn t < 


np 


57JO f 0* ) 9S1 


.0077 


np 


U4$jy/ 


.1} 1 l (.04-49) 


3 X 10." 7 


0009 


fi5<\r» r 01761 


5 * 10"* 


.0002 . 




.4:>44 (.034S) 


.9013 


n P 


/DO ^.U_U I ) 


.SI 20 


np 


D4S2965 


.5296 (.05S1) 


.306$ 


np 


/ niAA\ 


.2340 


np 


D4SI92 


.3 13d (.0406) 


.3715 


np 




.4525 


np 


U 4^420 


.:o9:> (,Ioj>V) 




np 




.1200 


np 


U451Q44 
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fi is the estimated proportion of alleles shared identical by descent, np, simulations not performed. 



individuals inherit disease susceptibility alleles but do not man- 
ifest the phenotype because of the simultaneous presence of 
protective alleles, ff model-based methods are used, it is impor- 
tant to provide a reasonably low estimate of penetrance and 
include a multifactorial component in the model. 

In the initial stages of analyzing a disorder like BPAD, which 
most likely displays multifactorial inheritance, robust model-free 
(allele sharing) methods are usually more useful than model- 
based linkage analysis (32). Concordant individuals should dem- 
onstrate excess allele sharing, even with the occurrence of phe- 
nocopies, genetic heterogeneity, high frequency of susceptibility 
alleles, and incomplete penetrance. Individuals who inherit sus- 
ceptibility alleles but do not manifest disease because of protec- 
tive alleles and individuals who inherit protective alleles but 
neve rt lie less manifest the disease will reduce the power of these 
analyses. Thus, regardless of the type of linkage analysis per- 
formed, the presence of protective alleles could have a major 
impact on identifying susceptibility loci. 

Although the idea that protective alleles could modify (or even 
prevent) a behavioral phenotype like BPAD is relatively novel, 
there are examples where such protective alleles can affect the 
expression or inheritance of other Mendelian and multifactorial 
disorders. The severity of sickle cell anemia is influenced by genes 
that increase the amount of circulating fetal hemoglobin (33). 
Similarly, the genotype of the chemokine receptor CCR5 dra- 
matically influences the kinetics of HI V-l infection, where mos4 
individuals who arc homozygous for a 32-bp deletion in the CCR5 
gene encoding the coreceptor for macrophage -tropic HI V-l are 
protected from virus infection (34). In Alzheimer's disease, 
apolipoprotein (Apo) E2, in contrast lo ApoE4, appears to 
reduce the risk of developing the disease and may protect 
individuals who inherit a disease-associated ApoE4 allele (35). In 
an extended Italian family, Apo A-Imilano protects against the 
development of both clinical and pathologic signs of atheroscle- 
rosis, despile significantly elevated plasma triglycerides and a 
markedly decreased level of high density lipoprotein-cholestcrol 
(36). In the nonobese diabetic mouse model of hum;wi autoim- 
mune insulin-dependent di:iheks mcllitus. pani;it protection 
from disease is pro-- iiled by resisianec allele* occurring >indy *'i 
either the KIJ3 or I.U1H) n..»:i-m:u>i hiM-«OMnp.it:!';!ii\ 



loci, whereas epistatic interaction between resistance alleles at 
these two loci produces nearly complete protection from diabetes 
(37). 

There are several mechanisms by which wellness or protective 
alleles could affect the clinical manifestations of BPAD in the Old 
Order Amish. One possibility is that dominant acting protective 
alleles, either singly or acting together in epistasis. could prevent 
or modify the BPAD phenotype. The variable penetrance of 
illness or its heterogeneous clinical manifestations could result 
from resistance or protective alleles that alone provide only 
partial protection, while together with other genes produce 
epistatic interactions, resulting in a grealer degree of modification 
of the phenotype. Alternatively, there also could be cellular target 
molecules, e.g., mood effectors, having forms that are either 
resistant or susceptible to the genetic and/or envtronmental 
susceptibility factors for BPAD. Individuals having resistant 
mood effectors would be protected from the effects of suscepti- 
bility alleles and/or environmental factors that result in the 
BPAD phenotype. In contrast, individuals with sensitive forms ol 
these mood effectors would be vulnerable to developing the BP1 
phenotype when requisite BPAD susceptibility alleles and/or 
environmental factors are present. If epistatic interactions arc- 
required for manifestation of the effects of either susceptibility or 
protective alleles, the existence of resistant and sensitive forms ol 
cellular effectors or protective alleles would be most apparent in 
families (or populations) where there is a high density of affected 
individuals such as the Old Order Amish in the present study. 
Regardless or the mechanism, the presence of wellness or pro- 
tective alleles would have a significant impact on linkage analyse* 
as evidenced by preventing the appearance of the BPAD phe- 
notype (or its presentation as a forme fruste) in individuals who 
are otherwise genetically predisposed to developing illness. 

Thus, we suggest that the identification of chromosomal loci 
harboring genes that contribute to the clinical manifestations ot 
BPAD will likely require a multilocus approach that considers 
both additive and subtractive influences of alleles on the Bl *u 
phenotype. The involvement ot protective or wellness alleles in 
'leterminin, the manifestation of the HP AD phenotvoc uouM 
|m ovule an aur:,eiKe. u-M..Hc CNnlanation f,n ;.t l-'"^ 1 " 1 ^; 1 
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Fig. 2. Model-free linkage analysis of wellness by using GH-PLUS. 
Map position is in Kosambi centi morgans. The -logio P was calculated 
by using P values gene rated by GH-PLUS on the assumption that the NPL 
score is standard normally distributed. A —login P of 4.0 corresponds 
asymptotically to a logarithm of odds score of 3.0. Only mentally healthy 
individuals 45 years of age or older were classified as being well (see text). 
(A) -logio P for markers on chromosome 4p: dottcdlinc, pedigree 110 
only: and solid line, pedigrees 110. 210. 310, and 410. (B) -tog l0 P for 
markers on chromosome 4q: doited line, pedigree 1 10 only, and solid line, 
pedigrees 110, 210, 310. and 4t0. 

The test statistics from our analyses for alleles linked to the 
absence of psychiatric illness in the Old Order Amish are at least 
as significant as those reported for any susceptibility locus, and 
further investigations into the significance of these findings in the 
inheritance of BPAD are warranted. The identification and 
characterization of protective alleles and their gene products 
could lead to the development of a more rational and direct 
approach to effective therapy for affective disorders. 
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The most characteristic features of bipolar affec- 
tive disorder (manic-depressive illness) are 
episodes of mania (bipolar I, BP1) or hypomania 
(bipolar II, BPII) interspersed with periods of 
depression. Manic-depressive illness afflicts 
about one percent of the population, and if 
untreated, is associated with an approximately 
20% risk of suicide 1 . Twin, family and adoption 
studies provide compelling evidence for a partial 
genetic aetiology, but the mode(s) of inheritance 
has not been identified 2 . Nonetheless, the majori- 
ty of genetic linkage studies have assumed classi- 
cal mendelian inheritance attributable to a single 
major gene. Although segregation analyses have 
yielded inconsistent results (with most studies 
rejecting a single locus inheritance model), the 
best single gene model is dominant inheritance if 
only BPI is considered 3 . Reported linkages of 
bipolar affective disorder on chromosomes 11, 
18, 21 and X have been difficult to substantiate, 
and additional studies are required for replication 
or exclusion of these regions 4 " 9 . We now present 
the results of our genome-wide linkage analyses 
that provide evidence that regions on chromo- 
somes 6, 13 and 15 harbour susceptibility loci for 
bipolar affective disorder, suggesting that bipolar 
affective disorder in the Old Order Amish is inher-" 
ited as a complex trait. 



The Amish study has employed rigorous standards 
for both clinical documentation and independent 
diagnostic evaluations with various reliability tests (see 
Methods, and Fig. 1). As the exact mode of inheritance 
ot bipolar affective disorder is unknown, we pei- 
formed linkage analyses using nonparametric methods 
(allele sharing, model-independent), as well as Iod- 
scorc analyses. Sixteen models were tested in the lod 
score analyses (see Methods). 

To evaluate the effect of analysing the data under 
multiple models, we compiled maximum lod scores 
for the 551 markers obtained under a single model 
(Fig. 2, Nl)(for dominant inheritance, nuclear fami- 
lies, homogeneity, and afTecteds only) and compared 
them with lod scores obtained under multiple testing 
(Fig. 2, N2). Mean lod scores were 0.18 versus 0.46, for 
Nl and N2, respectively. Thus, multiple testing raised 
lod scores on the average by 0.28 units. Also, no lod 
score exceeded 2.0 under the fixed model (Fig. 2, Nl), 
whereas three lod scores had values of at least 2.0 
under multiple testing (Fig. 2, N2). For the sib-pair 
analyses (see below), we used two weighting schemes 
and examined four categories of affection status. Also, 
the SIBPAL program furnished analyses for the mean 
number of alleles sharing identity by descent (IBD) in 
affected sib pairs as well as for regression of the num- 
ber of alleles sharing IBD on the squared difference of 
phenotypes in sib pairs, where phenotype = 0 for 
unaffected and phenotype = I for affected. This degree 
of multiple testing (12 tests) increases the type 1 error. 
As tests are non-independent it is difficult to judge the 
inflationary effect of multiple testing but the effect is 
likely to be moderate as shown by the more formal 
analysis of multiple lod score tests given above. 

We also performed analyses using the affected pedi- 
gree member (APM) method and the transmission 
disequilibrium test (TDT). The APM method fur- 
nished a few isolated low empirical P values that were 
not supported by other analyses. Thus, no APM results 
are given here. To assess a possible role of linkage dise- 
quilibrium, P values for the TDT were determined. 

Eleven percent of the markers (62/551) used in our 
genome- wide search gave a maximum lod score of 1.0 
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Fig. 1 Ancestral trace for Amish study bipolar 
pedigrees, a, Summary of each pedigree in 
'block' form, illustrating that all of the 8PI 
pedigrees trace along pathways leading to a 
common progenitor, one of some 30 couples 
that founded the present Lancaster County, 
Old Order Amish group 36 . This ancestral line 
encompasses the earliest cases of recurrent, 
psychiatric illness and the first confirmed 
cases of bipolar affective disorder. Bipolar I 
disorder among descendants of other pio- 
neers usually occurred after intermarriage 
with the BP progenitor line 14 . Note the LEFT 
extension coupled with the CORE Pedigree 
1 10, which provided the resource first used to report genetic linkage 
data 4 . After follow-up and addition of a RIGHT extension to Pedigree 
1 10, further genetic analyses were reported in 1989 (rat 5). Pedigree 
210 and partial Pedigree 310 (NIGMS Family 1075) 19 became a sec- 
ond large lateral extension 37 . Our report utilizes all of these earlier 
subjects plus additional expansions, especially in Pedigree 310. b, 
Summary of the diagnoses for the individuals in this study. One unaf- 
fected subject is counted twice — as a sibling/aunt of bipolar affect- 
eds in Pedigree 210, and as a spouse/mother of bipolar children in 
Pedigree 3 10. These Old Order Amish kinships continue to provide 
for lateral and lineal expansion and have evolved into the IX-Xth gen- 
erations of descendants at risk. 
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CLASS BOUNDS 
MAX LOD SCORES 



F<g. 2 Maximum lod scores, using 8PI as affected, for two scenarios: <N1) nuclear families homogeneity 
dominant inheritance, and affecteds only; and (N2) 16 combinations of analyses for each marker includinq 
dominant vs recessive inheritance, 5 pedigrees vs. nuclear families, homogeneity vs. heterogeneity and 
affecteds only vs affecteds and unaffecteds. N1 and N2 represent the number of markers furnishing maximum 
tod scores within given class boundaries for scenarios 1 and 2, respectively. 



or higher in at least one of the 16 analysis models. 
Consequently, only regions which yielded stronger evi- 
dence of linkage were considered further, namely those 
with markers having a test statistic of P<0.01 in any 
one analysis type, or maximum lod score of at least 
2.0. Six markers (D1S48. D6S7, D7S67, DUS146, 
DISS I and D15S45) met these criteria. Markers at 
three loci gave P values of 0.001 or smaller under SIB- 
PAL analyses: D6S7 (chromosome 6pter-p24) with P= 
0.0003; D13S1 {chromosome 13ql3) with P- 0.0003; 
and D15S45 (chromosome 15qll-qter) with P = 

0. 0003. Test statistics for these three markers, as well as 
for markers flanking these regions, are shown in Table 

1. Some of the formal P values furnished by the SIB- 
PAL program were re-estimated using computer simu- 
lations (3,000 replicates) run under an assumption of 
no linkage and are indicated by P 3 (see Methods). 
Therefore, they are not flawed by analyses that would 
furnish spuriously small formal P values. In lod-score 
analysis none of the markers reached lod = 3 criterion. 
However, in some of the nonparametric analysis meth- 
ods, P values of less than 0.001 (and even 0.0003, 
which is asymptotically equivalent to Z max = 2.5) were 
found. 

Results are typically stronger for BPI than for more 
liberal diagnostic categories (Table 1); that is, extend- 
ing the pool of affected individuals to include addi- 
tional psychiatric illness (BPII and recurrent MDD) 
appears to add *noisc* to the analyses. Generally, we 
obtain equivalent results for lod-score analyses and 
our ASP analysis, whereas the Haseman-Elston 
approach (SIBPAL program) typically provided 
stronger results. The main difference between the two 
programs, SIBPAL and ASP, is thai iherc is no weight- 
ing of multiple sib pairs in a sibship in SIBPAL. 

The relationship between point wise (locus -specific 



or nominal) and genome-wide significance levels was 
recently discussed by Lander and Kruglyak 10 . For sib- 
pair methods, pointwise P values 10 of 0.00074 and 
0.000022 correspond to suggestive and significant 
linkage, respectively, with 'significant* denoting a 
genome-wide P value of 0.05. For lod score analysis, 
the respective lod-score thresholds are 1.9 and 3.3. 

For D6S7, D13S1 and D15S45 (see Table I) an analy- 
sis of allele frequencies was carried out on probands 
and their mates to test for the existence of linkage dise- 
quilibrium 11 . Intraclass correlations for alleles within 
individuals and between individuals (who are mates) 
were calculated. In no case was there a significant dif- 
ference at the 5% level. For these three loci, as well as 
for the two markers flanking each of them, we tested 
for the presence of specific haplotypes. Again, no hap- 
lotypes were significantly more frequent than expected 
by chance when tested at the 5% level. 

A recendy described approach 12 for estimating the 
proportion of alleles sharing IBD based on all markers 
of a chromosome, was applied to the chromosomes 
carrying our best three markers (Table 2). For allele 
sharing among all possible affected sib pairs, it resulted 
in suggestive linkage for locus D6S7. However, TDT 
did not provide P values suggestive of linkage when 
applied to our best three loci and BPI affecteds (Table 

According to recently proposed criteria 10 , D6S7, 
D13S1 and D15S45 yield locus-specific P values sug- 
gestive of linkage. Therefore, it is imperative that the 
evidence for linkage of bipolar affective disorder to 
DNA markers on chromosomes 6, 13 and 15 in the 
Old Order Amish observed here be confirmed and 
extended in other pedigrees. Wc are collecting addi- 
tional members of the Amish pedigrees for analysis of 
more markers in the critical chromosomal regions. 
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Table 1 Results of linkage analyses 



Map 
distance (cM) 



Z^.dom Zro.rec SI8PAL P ASPP 



Table 2 Test results for TDT and multipoint affected ~ 



BPI 

D6S344 
D6S7 
D6SS9 
06S23 

8P1+ll b 
D6S344 
D6S7 
D6S89 
D6S28 

BP1+II 
D6S344 
D6S7 
D6SS9 
D6S28 

BPt+ll+MDO 
D6S344 
06S7 
D6S89 
D6S28 



0.0 
36.0 
17.7 



0.0 
36.0 
17.7 



0.0 
36.0 
17.7 



0.0 
36.0 
17.7 



.000 
2,342 
.097 
.000 

.000 
2.469 
.167 
.000 

.000 
1.885 

.732 
.000 

.000 
1.606 
.732 
.000 



.000 
1.456 
.001 
.003 

.000 
1.609 

.000 
.000 

.000 
.984 

.394 
-000 

.000 
.795 
.399 
.000 



.8729 
.0003* 
.5567 
.6262 

.9249 
.0000 
.7113 
.6272 

.9249 
.0003* 

.5241 
.6874 

.8977 
.0001 
.7126 
.6877 



BPJ 

D13S221 

D13S171 *l% 

D13S1 J'f 

D13S218 5- 

D13S291 101 

Bpi+n b 

D13S221 

D13S171 ™i 

D13S1 " 
D13S218 

D13S291 t0J 

BPI+ll 

D13S221 ... 

D13S171 } li 

D13S1 f 2 

D13S218 5- 

D13S291 iai 

BPkll+MDO 

D13S221 4 _- 

D13S1 f? 

013S218 5- J 

D13S291 101 



.2145 
.0513 
1.0000 
.7800 

.2341 
.0293 
.3230 
.7867 

.2036 
.1561 
.2743 
.6892 

.1055 
.4195 
.1714 
.7918 



.000 
.000 
1.276 
.312 
.056 

.000 
.000 
1-402 
.494 
.004 

.000 
.000 
1-203 

.307 
.006 

.000 
.000 
.000 
.000 
.000 



.012 
.102 

1.24a 

.664 
.175 

.000 
.000 
1.036 
.423 
.178 

.000 
.000 
.676 
.314 
.194 

.000 
.204 
-043 
.025 
.008 



1.0000 
.4905 
.0003* 
.0171 
.0865 

1.0000 
1.0000 
jOOOO 
.0175 
.1475 

1.0000 
1.0000 
.0090* 
.1403 
.2384 

1.0000 
1.0000 
.0223 
.3007 
.2262 



.9725 
.3736 
.0057 
.0641 
.3028 

.8902 
.4876 
.0056 
.0766 
.3998 

.7344 
.4665 
.0162 
.1484 
.4033 

.7672 
.9429 
.0932 
.3202 
.3525 





TDT 

(P value) 


Sib-pafr analysis (lod score) 
Indep. pairs A IJ pairs 


D6S344 
D6S7 
D6S89 
D6S28 


0.063 
0.200 
c\ *?nn 

0.500 


0.089 
0.920 
0 
0 


2.139 
2.211 
0.091 
0.039 


D13S221 

D13S171 

D13S1 

D13S218 

D13S291 


0.500 
0.500 
0.500 
0.497 

o!soo 


0 

0.598 
0.868 
0.595 
0 


0.002 
0.772 
1.422 
1.490 
0.375 


D15S49 

D15S48 

D15S45 

D15S117 

D15S148 


0.500 
0.500 
0.166 
0.500 
0.119 


0.737 
0.239 
0.877 
0.084 
0.073 


0.517 
0.247 
1.540 
0,661 
0.629 



BPI 

D15S45 , c 
D15S117 J-J 
D15S148 
D15S38 
D1SS36 

D1SS4S 

D15S117 J'S 

D15S148 

D15S38 

D15S36 
BPWI 

D1SS4S K . 

D15S117 *S 

D15SU8 

D15S38 

D1SS38 
BPI+lkMOO 

D15S4S 

D15S117 

D15SU8 

D1SS38 

D15S36 



1.2 
6.1 
0.0 



1.2 
6.1 
0.0 



1.2 
6.1 
0.0 



5.6 
1.2 
6.1 
0.0 



1.114 

.130 
.338 
.000 
.355 

1.097 
.332 
.752 
.067 
.646 

.857 

.089 
.461 
.000 
.368 

1.709 

.000 
.148 
.000 
.000 



.798 
.580 
.610 
.000 
.400 

.446 

.589 
.613 
.000 
.402 

.731 
.726 
.829 
.000 
.292 

.473 

.096 
.192 
.000 
.000 



.0003 a 

.0843 
.0217 
1.0000 
.0114 

.0018 
.1225 
.0219 
1.0000 
.0118 

.0183 

.0910 
.0123 
1.0000 
.0131 

.0032 

.2119 
.0360 
1.0000 
.1423 



.0163 

.1660 
.1000 
.9853 
.0662 

.0456 
.2346 
.0976 
.9904 
.0844 

.0399 

.1825 
.0589 
.8551 
.1172 

0)150 
.3546 
.0998 
.7914 
■3168 



Zm^,: maximum tod score in analysis o) nuclear families, affected only, pene- 
trance ratio (genetic versus nongenetic cases) of 500:1. with allowance (or het- 
erogeneity (exception: for 06S7, afJocteds and unaffecteds) 
f^T 0 ™ ^ **** roc: under dominant or recessive Inheritance 
S1BPAL P: P values furnished by SfBPAL program in t-test for excess allele shar- 
IntMn affected sib pairs (exception: resufta for regression analysis given for 

ASP p\ P values ! in t-test for excess of allele sharing in aflected sibs, multiple s.b 
pairs in same slbshlp weighed by number of metoses 
Clmlcal categories: MOO includes onry recurrent major depressive disorder 
Number of affected* in clinical hiorachies were: 31 BPI. 35 BPUBPtl 1 ' 39 ' 
BPI+BPII, and 49 BPI+BPII+MDO 

•Empirical P vahjes were estimated by computer simulation 

. 6P " 00X9 borderline BP I are Included (such as clinical BP 

I and ROC manic) 

Past attempts to replicate linkage findings for bipo- 
lar affective disorder have been plagued by diagnostic 
uncertainties, genetic heterogeneity, phenocopies, 
genotyping errors and the complexities of performing 
and interpreting statistical results. Our study, however, 
represents a 19-year longitudinal study of an isolated 
population in which there is a relatively narrow spec- 



MAKER/siBS program - pairing options 2 (independent pairs) and 3 
(all possible pairs), with allowance for dominance variance. 

trum of illness (not one case of schizophrenia occurs 
in these . pedigrees); bipolar disorder is the predomi- 
nant diagnosis. The rigorous longitudinal assessment 
of these pedigrees combined with the systematic and 
blind psychiatric evaluations and diagnoses should 
also greatly reduce the number of misdiagnoses 
included in our linkage analyses. Moreover, the 
restricted gene pool characteristic of this relatively 
closed population should reduce the number of dis- 
ease-causing alleles, minimizing the problem of genet- 
ic heterogeneity. 

Similar to other common and complex diseases like 
diabetes, hypertension, and perhaps even schizophre- 
nia, our data suggest that genetic factors likely con- 
tribute to the pathogenesis of bipolar affective 
disorder, where in the majority of cases, inheritance is 
multifactorial rather than simple mendelian transmis- 
sion. Like the genetic variance observed for quantita- 
tive traits, bipolar affective disorder (even in a relative 
genetic isolate like the Old Order Amish) appears to be 
a polygenic (complex) trait resulting from the variable 
effects of multiple genes. The results of our genome 
wide scan suggest that genes on chromosomes 6, 13 
and 15, rather than just different mutant alleles of a 
single gene, determine the susceptibility to and pheno- 
type of bipolar affective disorder in the Old Order 
Amish. Additional sets of genes may underlie the sus- 
ceptibility to develop bipolar affective disorder in 
other populations. Although we emphasize the provi- 
sional nature of these loci, the findings from our 
genome wide search in the Amish can now be exam- 
ined in other populations, and may eventually lead to 
the identification of specific genes responsible for 
bipolar affective disorder. 

Methods 

Diagnostic assessment. The genetic -epidemiologic study of 
bipolar affective disorder among (he Old Order Amish in 
southeastern Pennsylvania has been described' 314 . Case ascer- 
lainmeni fur mental illness among the Amish began with a 
community-wide network of informants and institutional ros- 
ters reviewed with informed consent 15 . Over 400 patient cases 
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have been ascertained. A psychiatric review board composed 
(since 1976) of lames N. Sussex, Abram M. Hostetter, John J. 
Schwab, David R. Offord and Jean Endicott used both psychi- 
atric interviews"' and abstracted medical records to perform 
diagnostic assessments based on strict research diagnostic cri- 
teria (RDC) 17 . Assessments by this review board were made 
blind to pedigree membership, diagnostic opinions and treat- 
ment information in the medical records, and genetic marker 
status. As the Board's diagnostic procedures yielded confirmed 
cases of BP! affective disorder, the immediate families of these 
patients were evaluated for psychopathology. Pedigree 1 10 was 
selected (in 1981) for initial genetic linkage study because of 
relationships between nuclear families, based on BPI probands, 
and illness spanning several generations'*. When one examines 
the relative risk for individuals used in our linkage study, there 
is a very high prevalence of affective disorder, with age-correct- 
ed morbid risk rates for BPI, BPII, and MDD (major depres- 
sive disorder) of 17%, 4% and 6%, respectively. This gives an 
overall rate of 27% for major affective disorder in our pedi- 
grees. The present sample, which includes extensions to the 
original family (Fig. \a) totals 207 members, with 31 diag- 
nosed BPI, 50 with other psychiatric diagnoses, and 126 unaf- 
fected individuals (Fig. \b). Over 125 medical records were 
abstracted and Board reviewed to document the 31 BPI cases. 
The average age of onset for BPI disorder was 22 years. Relia- 
bility of the bipolar diagnoses was checked when 16 of the 31 
cases (52%) were evaluated twice, with an average five year 
interval between the blind assessments using different clinical 
documentation and resulting in 100% concordance. The high 
reliability obtained lessens the likelihood of misdiagnoses or a 
false positive BPIin our linkage analyses 18 . 

Apart from RDC diagnoses, the project psychiatric panel 
also recorded clinical opinions in a consensus 'clinical diagno- 
sis'. There was 100% concordance between these two types of 
diagnostic conclusions (5 board members) for the 31 BPI cases 
and 13 of the 15 cases of recurrent major depressive disorder 
(Fig. \b). Of particular interest are the diagnostic results for 
the 8 cases of BPII. Four were designated BPII by both RDC 
and clinical opinion. The other four were labelled BPI accord- 
ing to clinical opinion, and two were classified as 'probable 
BPI' by the strict RDC This is important because it shows that 
true BPII disorder occurs rarely in these pedigrees; BPII 
appears more as a *BPP waiting to happen. 

Patient samples. Only after the initial diagnostic assessments of 
probands and their relatives were the blood samples obtained 
from the subjects using informed consent Codes were used to 
maintain confidentiality. The phlebotomist was kept blind to 
pedigree relationships and diagnostic status. Given a variable 
age of onset for the affective disorders and the impact of new 
diagnoses on data analyses, the Amish Study policy is to con- 
duct annual follow-up. This permits an update of any changes 
in the course of illness for affected individuals and allows us to 
see if previously well individuals have onset with illness. It also 
facilities ascertainment of new cases of major affective disorder 
in close branches of the family so plans can be made to include 
them in pedigree expansion. Blood samples were collected with 
informed consent and lymphoblastoid cell lines were estab- 
lished at the Coriell Institute of Medical Research and/or the 
National Institute of Mental Health. The NIGMS Human 
Genetic Mutant Cell Repository catalog contains updated pedi- 
gree and diagnostic information 19 . 

Genotyping. DNA was extracted from peripheral blood sam- 
ples and/or immortalized lymphoblastoid cell lines 20 . Mapping 
panels were constructed to determine the best order of markers 
typed on the bipolar pedigrees using genotypic data from the 
CEPH version 7 database, using the MultiMap linkage analysis 
program 21 . The RFLP and microsatellilc markers used resulted 
in a linkage map with an average spacing of between 5 and 10 
cM 222 \ Microsalcllitc markers were genotyped individually as 



described 9 -*, and by an adapted multiplex procedure'* In Hie 
multiplex genotyping procedures, a toud of 25 markers were 
analysed in each lane of the genotyping g c!s . To accoim ,| ish 
this, five markers were coamphfied in each PCR lube and five 
sets of five markers were pooled and precipitated prior to gel 
loading. Sets of five markers were amplified in 20 pi rcaclions 
that included: 10 mM Tris-HCl, pH 8.3, 50 mM KC1 0 001% 
gelatin, 1 .5 mM MgCl 2( 0.2 mM dNTPS. 0.5 U Taq polymerase. 
I pM of each primer (10 primers in total), and 50 ng genomic 
DNA. Samples were denatured at 94 °C for 1 min, followed by 
30-35 cycles at 94 °C for 1 5 s, 55 °C for 1 5 s, and 72 *C for 1 5 s. 
After the final cycle, the reactions were incubated at 72 °C for 3 
min. Following amplification, five sets of PCR amplifications 
were pooled and isopropanol precipitated in 96 well microtiter 
plates. The pooled PCR products were rcsuspended in 10 pi 
loading dye containing formamide, denatured at 95 °C for 5 
min -and loaded onto 6.0% denaturing polyacrylamide gels. 
After electrophoresis, the gel fractionated microsatellilc mark- 
ers were transferred to nylon membranes by capillary transfer, 
and visualized by hybridizing membranes with marker specific, 
chemiluminescent probes. To make the probes, one of the 
oligonucleotides used to amplify each marker was labelled 
using the ECL detection kit (Amersham). Multiple probes cor- 
responding to markers of different sizes were hybridized to the 
filters simultaneously. Chemiluminscent signals were detected 
by autoradiography. Allele sizes for the microsatellite markers 
were determined relative to a pUC18 sequence or Sequamark 
marker ladder (Research Genetics). To maintain allelic desig- 
nations for the purpose of allele frequency calculations, DNA 
samples from replicate individuals were included within and 
between gels. Films were scored either manually or using semi- 
automated atlele-calling software (Bio Image), and were inde- 
pendently analysed by two individuals blind to disease status. 
Data from all markers were transferred into the same file sys- 
tem for linkage analyses. 

linkage analyses. Statistical analyses for linkage were per- 
formed using both nonparametric (SAGE SlBPAL program for 
Haseman-Elston sib pair tests; affected sibpair analysis (ASP) 
with weighting of multiple affected in the same subship by 
number of meioses; affected pedigree member (APM); and the 
transmission disequilibrium test (TDT)) as well as lod score 
analyses. Sixteen models were tested in the lod score analyses, 
and included: dominant versus recessive inheritance; large 
pedigrees versus data broken into nuclear families; homogene- 
ity versus allowing for heterogeneity; and analysis for affectcds 
only versus using affecteds and una/fecteds. For nonparametric 
analyses, based on the asymptotic (theoretical) distribution of 
the test statistic the SIBPAL program 26 furnished formal P val- 
ues in the test for an excess proportion of alleles shared IBD 
and in the Haseman-Elston regression test 27 . The 'true P value 
(the empirical significance level) is defined as the probability 
that the observed result, or more extreme one, is obtained by 
chance alone. To estimate empirical significance levels associat- 
ed with these results rather than relying on die formal P value, 
we carried out computer simulations (3,000 replicates) for 
each marker with a formal P value of 0.01 or less. These simu- 
lations were extremely time consuming as a complete analysis 
had to be carried out for each replicate. Therefore, empirical P 
values were obtained only for those analyses on the BPI and 
BPI + II clinical diagnostic categories where markers had a for- 
mal P value of 0.001 or less. Resulting empirical P values (iden- 
tified by P 1 in Table 1) typically were about ten times higher 
than the formal 'P values' issued by the program. In addition, 
sib-pair analyses were performed using a computer program in 
which multiple pairs per sibship are weighted by the number of 
meioscs, and that also carried out the TDT analyses 28-31 . Mul- 
tiple alleles were taken into account as described in an exten- 
sion to the TDT- 12 . The TDT was applied to allow for a possible 
involvement of linkage disequilibrium. The total number of 
affected sib pairs in our data set was 12 for BPI, 13 fur BPI + 
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BPII 1 (borderline BPI), 19 for BPI + BPI1, and 30 for BPI + 
BPH + MDD. While those numbers arc small, they do not rep- 
resent random sib pairs but rather come from n unique closed 
population and are, presumably, all at least distantly related to 
each other. 

Allele sharing based on all markers on a chromosome was 
applied as implemented in the MAPMAKER/S1DS program, 
resulting in lod scores along points on the chromosome. Two 
program options for forming affected sib pairs were used, 
'independent pairs' (first sib versus each of the other sibs) and 
'all pairs'. However, because sib pairs arc unweighted, the latter 
option may furnish lod scores that are too high in small to 
moderate sample sizes. 

For parametric analyses, two -point lod scores were calculat- 
ed with the LINKAGE programs 3 *. Each marker was analysed 
under 16 models (dominant versus recessive inheritance, large 
pedigrees versus data broken into nuclear families, homogene- 
ity versus allowing for heterogeneity by the HOMOG pro- 
gram 3 * 1 , and for afTecteds only versus affected and unaffected 
individuals considered). A penetrance ratio (genetic versus 
non-genetic cases) of 500:1 and disease allele frequencies 
adjusted to reflect a population prevalence of 1% were used. 
Individuals without psychiatric illness under any diagnostic 
scheme were considered unaffected, whereas those not catego- 
rized as affected under one scheme but affected under one of 
the other diagnostic categories were taken to be unknown. 



Linkage disequilibrium influences certain types of identity- 

t y ih Zl i cr h sib " p3ir strat T * The basi ™» ^ 

lh * aIldc sharm 5 reIaliv <* in isolated populations 

may be exaggerated due to the fact that large regions of 
founder chromosomes have not been broken up by rcconibina 
non. For this reason, we employed the TDT to test for link-ice 
and applied IBD methods, which are unlikely to be significant 
influenced by disequilibrium. 
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Melanocortin-4 receptor gene variant il03 is negatively associated with 
obesity . 

Geller F, Reichwald K, Dempfle A, Illig T, Vollmert C, Herpertz S, 
Siffert W, 

Platzer M, Hess C, Gudermann T, Biebermann H, Wichmann HE, Schafer H, 
Hinney A, 
Hebebrand J. 

Institute of Medical Biometry and Epidemiology, Philipps -University , 
Marburg , 
Germany . 

Several rare mutations in the melanocortin-4 receptor gene (MC4R) 
predispose to 

obesity. For the most common missense variant V103I (rs2229616) , 
however, the 

previously reported similar carrier frequencies in obese and nonobese 
individuals are in line with in vitro studies, which have not shown a 
functional 

implication of this variant. In the present study, we initially 
performed a 

transmission/disequilibrium test on 52 0 trios with obesity, and we 
observed a 

lower transmission rate of the 1103 allele (P=.017), which was an 
unexpected 

finding. Therefore, we initiated two large case-control studies 
(N=2,334 and 

N=661) and combined the data with those from 12 published studies, for 
a total 

of 7,713 individuals. The resulting meta-analysis provides evidence for 
a 

negative association of the 1103 allele with obesity (odds ratio 0.69; 
95% 

confidence interval 0.50-0.96; P=.03), mainly comprising samples of 
European 

origin. Additional screening of four other ethnic groups showed 
comparable 1103 

carrier frequencies well below 10%. Genomic sequencing of the MC4R gene 
revealed 

three polymorphisms in the noncoding region that displayed strong 
linkage 

disequilibrium with V103I. In our functional in vitro assays, the 
variant was 

indistinguishable from the wild- type allele, as was the result in 
previous 

studies. This report on an SNP/haplotype that is negatively associated 
with 

obesity expands the successful application of meta-analysis of modest 
effects in 

common diseases to a variant with a carrier frequency well below 10%. 
The 

respective protective effect against obesity implies that variation in 
the MC4R 

gene entails both loss and gain of function. 
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Apolipoprotein-E s4 Allele Frequency and 
Conferred Risk for Cuban Americans With 

Alzheimer's Disease 

Steven Sevush, , M.D., Gloria Peruyera, , B.A. 
Fiona Crawford, , Ph.D., Michael Mullan, , M.D., Ph.D. 



The authors compared Apo-E e4 frequencies for Cuban American and white non- 
Hispanic persons with and without Alzheimer's disease (AD). The Apo-E e4 allele con- 
ferred as large a risk for AD in Cuban Americans as for white non-Hispanics, a finding 
that differs from those for Hispanic subjects as a whole. (Am J Geriatr Psychiatry 2000; 
8:254-256) 



An association between Alzheimer's disease (AD) 
and the apolipoprotein-E (Apo-E) e4 allele has been 
established worldwide. 1 The strength of the associa- 
tion, usually measured by comparing Apo-E e4 allele fre- 
quencies for patients with AD with those for age- 
matched control subjects without dementia, has been 
reported to vary across ethnicities and geographical ar- 
eas. 2 The ethnic variations in Apo-E e4-conferred risk 
for AD are potentially important because they may cor- 
respond to ethnic differences in pathogenetic mecha- 
nisms 1 or in clinical manifestations of the disease. 3 

Cross-ethnic comparisons of Apo-E e4 allele fre- 
quency in patients with AD in the United States have 
typically utilized a three-way ethnic partition, distin- 
guishing white non-Hispanics, Hispanics, and African 
Americans as distinct subgroups. Using this scheme, as- 
sessments of Apo-E e^conferred risk for AD by Tang et 
al. 4 and Farrer et al. 2 revealed a lower risk for Hispanic 
than for white non-Hispanic individuals. More recently, 
Tang et al., 1 utilizing data from a prospective longitu- 
dinal study, made the stronger argument that the Apo- 



E e4 allele conferred no risk of AD at all for Hispanic 
individuals. The low or absent Apo-E e4-associated risk 
for AD reported for Hispanic individuals may not, how- 
ever, be generalizable to all Hispanics but may be dif- 
ferent for different Hispanic subgroups. 

In the present study, we took advantage of the avail- 
ability of a homogeneous Hispanic subgroup compris- 
ing Cuban American individuals living in South Florida. 
We evaluated Apo-E genotypes for Cuban American and 
white non-Hispanic patients with AD and age-matched 
control subjects without dementia and compared the 
risk for AD conferred by the Apo-E e4 allele for both 
groups. 



METHODS 

A total of 145 patients with AD (80 Cuban American 
and 65 white non-Hispanic) and 49 control subjects (21 
Cuban American and 28 white non-Hispanic) served as 
the study group. Patients with AD were chosen consec- 
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utively from patients diagnosed with probable AD ac- 
cording to NINCDS criteria 5 after a comprehensive clini- 
cal and laboratory assessment conducted by the State 
of Florida-funded University of Miami Memory Disor- 
ders Clinic. Control subjects were either spouses (83%) 
or friends (17%) of patients. Patient ethnicity was de- 
termined by self-report (or caregiver report) of ethnicity 
("Hispanic" or "white non-Hispanic") and country of 
birth. All subjects in the Cuban American group classi- 
fied themselves as Hispanic and gave their country of 
birth as Cuba. Age at disease onset for AD patients was 
estimated via structured interview of the patient's care- 
giver by two independent examiners (interrater Pearson 
r = 0.85). The average of the two estimates was used for 
analyses. 

Apo-E genotypes were determined by PCR amplifi- 
cation according to the method previously described by 
Bennett et al. 6 Allele frequencies were determined by 
counting alleles and calculating proportions relative to 
the total number of alleles in the sample. Chi-square was 
used to compare Apo-E e4 allele frequencies for patients 
with AD with those for control subjects for both the 
Cuban American and white non-Hispanic groups. To de- 
termine whether the increase in Apo-E e4 frequency for 
patients with AD relative to controls was different be- 
tween ethnicities, a method for comparing odds ratios 
described by Cohen 7 was used. First, phi-coefficients 
(<J>; i.e., the fourfold point correlation coefficient) were 
calculated for the contingency tables corresponding to 
each ethnicity; then a Fisher z transformation was ap- 
plied; and finally, the statistical significance of the dif- 
ference between the transformed correlation coeffi- 



TABLE 1. Apo-E genotype and allele frequencies, n (percent) 

Cuban American 



Genotype 






4/4 


6(7.5) 


0(0) 


4/3 


28 (35.0) 


2 (9.5) 


4/2 


0(0) 


1 (4.8) 


3/3 


40 (50.0) 


17(81.0) 


3/2 


7(7.5) 


1 (4.7) 


2/2 


0(0) 


0(0) 


Allele frequency 






2 


0.037 


0.48 


3 


0.713 


0.881 


4 


0.250 


0.071 


Odds ratio 




4.33 


95% Confidence limits 




1.27-14.79 


Note: AD = Alzheimer's disease. 





cients was evaluated with TV set as the harmonic mean 
of the ns for the two groups. 



RESULTS 

There were no significant differences in age or sex be- 
tween Cuban American and white non-Hispanic sub- 
jects for either patients or control subjects. Also, there 
was no significant difference between Cuban American 
and white non-Hispanic patients with AD for estimated 
age at disease onset. Patients and control subjects were 
matched for age, but patients with AD were more likely 
than control subjects to be female (P- 0.005). 

The Apo-E genotypes and allele frequencies for Cu- 
ban American and white non-Hispanic patients with AD 
and control subjects are shown in Table 1. Apo-E e4 
allele frequencies were greater for patients with AD 
than for control subjects, both for white non-Hispanic 
(P = 0.005) and for Cuban American subjects (P= 0.01). 
The relative risk of AD conferred by the Apo-E e4 allele 
was somewhat greater for Cuban Americans (OR = 4.33; 
range: 1.27-14.79) than for white non-Hispanics 
(OR = 3.40; range: 1.48-7.79), but the difference was 
not statistically significant (P>0.1). 



DISCUSSION 

For both Cuban American and white non-Hispanic 
groups, the Apo-E e4 allele frequency was greater for 



White Non-Hispanic 



8 (12.3) 


1 (3.6) 


28 (43.1) 


5 (17.9) 


3 (4.6) 


1 (3.6) 


19(29.2) 


18 (64.3) 


7 (10.8) 


3 (10.6) 


0(0) 


0(0) 


0.077 


0.071 


0.562 


0.786 


0.361 


0.143 



3.40 
1.48-7.79 



AD (n = 80) Control Subjects (n = 21) AD (n = 65) Control Subjects (« = 28) 
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patients with AD than for control subjects without de- 
mentia, indicating that the Apo-E e4 allele conferred a 
risk for AD for both groups. Comparison of the phi- 
coefficients for the corresponding contingency tables 
indicated that the risk was not significantly different be- 
tween the two groups. The substantial risk of AD con- 
ferred by the Apo-E e4 allele for Cuban Americans in 
this study contrasts with the previously reported low or 
nonexistent risk of AD associated with the Apo-E e4 
allele for Hispanics in North Manhattan 1,4 or of unspe- 
cified composition. 2 It is possible that differences in the 
ethnic makeup of the Hispanic groups in those studies 
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and the present one are at least partly responsible for 
the differences in findings. Schur et al. 8 have noted that 
it is important to distinguish among such Hispanic sub- 
groups as Cuban American, Mexican American, and His- 
panics of Puerto Rican or Dominican origin when con- 
ducting clinical studies comparing Hispanic and 
non-Hispanic groups. The findings in the present study 
suggest that such a distinction may be necessary when 
assessing the implications of Apo-E e4 allele frequency 
and conferred risk for AD. 
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Alzheimer disease (AD) risk is significantly influenced by the 
APOE2 and APOE4 alleles. In turn, the -49IATand THI/E47cs poly- 
morphisms alter APOE gene expression levels. To determine 
whether these two alleles exert any significant effect on AD devel- 
opment we have analysed the genotypes of the APOE promoter 
-491 AT and Thl/E47cs polymorphisms in 163 AD patients and 155 
controls divided into three age at onset/age dependent subgroups. 



Our study has detected aThl/E47cs-Tallele accumulation in healthy 
individuals over 75 years of age, which suggests it plays a protective 
role against AD.TheThl/E47cs-Tallele may provide greater protec- 
tion against AD than APOE2, although this awaits proof of Thl/ 
E47cs-T allele overrepresentation in healthy individuals of other 
populations. NeuroReport 13:1403-1405 © 2002 Lippincott 
Williams & Wilkins. 
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INTRODUCTION 

The apolipoprotein E (APOE) gene is known to be a major 
risk factor for AD in many populations [1,2]: the APOE4 
allele increases the risk of AD in a dose-dependent manner 
[3,4] and the APOE4/4 genotype may precipitate the onset 
of the disease [4]. On the other hand, since the APOE2 allele 
shows a high prevalence among centenarians [5,6], some 
authors have concluded that it protects against AD [7,8]. 
Statistically, whereas AD will develop in 50% of APOE4/ 4 
homozygotes before they reach 70 years of age, the median 
age at onset of AD among APOE2/3 heterozygotes is over 
90 years [9]. Moreover, in contrast to the fast rate of 
cognitive decline shown by APOE4/4 patients, the APOE2 
allele is associated with slow AD progression [10]. 

A flanking markers study has suggested that polymorph- 
isms in the APOE locus itself may also confer risk for 
developing AD [11]. This hypothesis has gained support 
from the observation of an allelic distortion in relative 
APOE mRNA levels in the brain of AD patients as 
compared with controls [12]. Several polymorphisms in 
the APOE promoter region have also been identified [13- 
15]. Two of these polymorphisms (-491 AT and Thl/E47cs) 
influence APOE gene expression levels, which are dimin- 
ished by the — 491-T allele and increased by the Thl/E47cs- 
G allele [14]. 

Different studies have reported a deleterious effect of the 
Thl /E47-T allele [15] and a protective effect of the -491-T 
allele, the -491 -AA genotype being responsible for in- 
creased risk in AD patients [13]. An additional study has 



confirmed the deleterious effect of the Thl/E47-T allele and 
the protective effect of the -491-T allele but only in APOE3/ 
4 genotype carriers [16]. If these alleles alter APOE gene 
transcription levels, the risk conferred by the APOEe4 allele 
or the protective effect exerted by the APOEe2 allele should 
vary accordingly [15]. 

Several authors have proposed the existence of protective 
genes against AD. In this regard, AD prevalence has been 
found to be decreased in a cohort with a high degree of 
Cherokee Indian ancestry studied in Oklahoma [17]. The 
low frequency of AD in this population was not changed by 
the APOEe4 allele. Similar data have been reported when 
the Cree Indian tribe was compared with an urban white 
population [18]. 

Since AD risk is significantly influenced by the APOE and 
APOE4 alleles and, in turn, APOE gene expression is altered 
by the polymorphism -491 AT and Thl/E47cs in the APOE 
gene promoter region, we have studied these two poly- 
morphisms in three age at onset/age dependent subgroups 
of AD patients with the aim of detecting differences in allele 
representation between patients and controls. 

MATERIALS AND METHODS 

The study included 163 AD patients (age range 48-85 years; 
mean age 70.8 years; male:female ratio, 1: 1.6) with clinical 
diagnosis of probable AD according to the DSM-IV and 
NINCDS-ADRDA criteria [19]. Also investigated were 155 
control subjects (age range 45-84 years; mean age 68.7 years; 
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male:female ratio, 1:1.5) with no neurodegenerative dis- 
order. Patients and controls were divided into subgroups 
according to their age at onset/age as follows: <65 years 
(46 patients, 46 controls), between 65 and 74 years (69 
patients, 65 controls), and > years (48 patients, 44 controls). 

DNA was extracted from peripheral blood cells according 
to standard procedures. The common APOE gene poly- 
morphism was studied as described previously [20]. A 
nested PCR was used for promoter polymorphism detec- 
tion. First, a fragment of 418 bp was amplified with primers 
APOE491A: CACCATgTTggCCAggCTggTCTCAA and 
TH1/E47L: ggAggAAggAggTggggCATAgAgg. PCR reac- 
tion mix contained 5mM MgCl 2 , and annealing and 
extension were carried out in a single step of lmin at 
70°C. The obtained fragment served as template for 
amplification of two fragments containing the -491AT and 
TM/E47 polymorphisms as described [14]. 

Allele frequencies were calculated by the allele counting 
method. The y 1 test of association was employed for the 
assessment of allele frequency differences between controls 
and patients, whereas calculation of odds ratios by logistic 
regression analysis was used to determine the effects of age, 
gender, and the -491 AT and Thl/E47cs polymorphisms on 
healthy controls. A and T allele frequency values were used 
for -491 AT polymorphism analysis, whereas the Thl/E47cs 
polymorphism was evaluated by means of T and G allele 
frequency values. Correction for the presence of both the 
-491-Tand Thl/E47cs-G alleles in APOE3/3 and APOE4/ 
4 homozygous individuals was applied. 



RESULTS 

When the AA, AA/AT, TT, or TT/AT genotype frequencies 
of the -491 AT polymorphism and the TT, TT/TG, GG, or 
GG/TG genotype frequencies of the Thl/E47cs polymorph- 
ism were analysed, the results obtained did not reveal any 
significant differences in genotype distribution between 
patients and controls for either the -491 polymorphism or 
Thl/E47 polymorphism. Taking into account that the -491- 
T allele diminishes APOE gene transcription and the Thl / 
E47-G allele enhances it, corrections were applied for the 
presence of both the -491-Tand Thl/E47-G alleles in APOE 
313 and APOE 4/4 homozygous carriers in our data 
analysis. 

The -491 AT polymorphism exhibited T-aliele frequen- 
cies lower than 0.1 in patients and between 0.09 and 0.13 in 
controls, with no significant differences among the various 
age at onset-dependent/ age subgroups (Table 1). 

TM/E47 polymorphism allele frequencies did not show 
significant differences between patients and controls in 
subgroups including individuals up to 74 years of age. On 
the contrary, the control subgroup formed by subjects of 75 
years or older exhibited a very high Thl/E47-T allele 
frequency when compared with patients and the other age- 
dependent control subgroups. Testing of the association 
between the TM/E47-T allele and lack of AD provided a 
high odds ratio for the subgroup consisting of individuals 
aged >75 years (Table 2). 

As expected, APOE allele frequencies did not present 
significant differences between cases and controls. Fre- 
quency-values were 0.02, 0.03 and 0.03 for the three AD 
subgroups and 0.06, 0.05, and 0.06 for the control subgroups. 



Table I. -49IAT polymorphism allele frequencies among AD cases and 
controls 



Group 


Number 


Alleles 










AD/C 


A (AD/C) 


T(AD/C) 




OR b 


<65 


46/46 


0.94/0.91 


0.06/0.09 


0.294 


n.s. 


65-74 


69/65 


0.95/0.90 


0.05/0.10 


0.096 


n.s. 


>75 


48/44 


0.93/0.87 


0.07/0.13 


0.124 


n.s. 



"Value for / 2 test of differences in -491 AT polymorphism allele frequencies 
between cases and controls. 

b Odds ratios (OR) estimates for the effect of the -491 -A allele on risk for 
AD, 



Table 2. Thl/E47cs polymorphism allele frequencies among AD cases and 
controls 



Group 


Number 


Alleles 










AD/C 


A (AD/C) 


T(AD/C) 




OR b 


<65 

65-74 

>75 


46/46 
69/65 
48/44 


0.59/0.55 
0.52/0.44 
0.48/0.83 


0.41/0.45 

0.48/0.56 

0.52/0.17 


0.421 
0.107 
< 0.001 


n.s. 
n.s. 

4.8(2.3-9.4) 



*Value for % 2 test of differences inThl/E47cs polymorphism allele frequen- 
cies between cases and controls. 

b Odds ratios (OR) estimates for the effect of theThl/E47-G allele on risk 
for AD. 



These values are in agreement with known APOE allele 
frequencies in AD patients and the normal population. 

DISCUSSION 

The results yielded by the present study allow us to propose 
the Thl/E47-T allele as a protective factor against AD in 
individuals >75 years of age. Our data are in disagreement 
with those obtained by Lambert et al., who have reported a 
deleterious effect for the same allele [15]. These discrepan- 
cies are probably due to differences in data analysis. 
Specifically, in the study by Lambert et al which included 
293 AD patients (age at onset 71.02 ±8.9 years; range 49-94) 
and 310 controls (age 73.5 + 10.9 years; range 59-105), 
division into age at onset /age dependent subgroups was 
not introduced and correction for the presence of both the 
-491-Tand TM/E47-G alleles in APOE3/3 and APOE4/4 
homozygous individuals was not applied (the existence of 
the -491 AT polymorphism was then unknown) [15]. 

Formation of amyloid plaques and neurofibrillary tangles, 
the neuropathological hallmarks of AD, precedes by many 
years the appearance of the first AD symptoms. In fact, 
these neuropathological changes have become irreversible 
by the time the disease is diagnosed [21,22]. Among the 
many molecules taking part in the formation of amyloid 
plaques, apoE outstands for being essential for p-amyloid 
deposition and insoluble fibril formation, as shown by data 
obtained in an APOE knockout mouse model studies 
[23,24]. Since APOE gene expression is inhibited by the 
Thl/E47-T allele and enhanced by the Thl/E47-G allele 
[14], apoE levels are lower in the brains of the rare Thl/E47- 
T allele carriers than in those of the commoner Thl /E47-G 
allele carriers. Decreased apoE levels would lead to 
diminished interaction between apoE and ^-amyloid in T 
allele carriers, whose resilience to fibril formation would 
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result in the slowing-down of amyloid plaque generation 
and, ultimately, increased resistance to AD development. 

Our present results need confirmation by studies of 
healthy subjects >90 years of age. If similar data are 
obtained in other populations, the Thl/E47cs-T allele would 
emerge as an important protective factor against AD. In fact, 
the protective role of the Thl/E47cs-T allele would be even 
greater than that of the APOE2 allele, since the former seems 
to exert its beneficial influence against AD development 
already at 75 years of age whereas the latter is unable to 
provide protection until individuals have surpassed 90 
years of age. 

CONCLUSIONS 

The Thl/E47cs-T allele of the APOE promoter region seems 
to be a protective factor against AD development in the 
Spanish population. This effect appears to be the result of 
lower APOE gene expression in Thl/E47cs-T allele carriers 
in comparison with the commoner Thl/E47cs-G allele 
carriers. Lower apoE accumulation rates in the brain would 
slow down amyloid plaque formation and, consequently, 
provide protection against AD development. 

The Thl/E47cs-T allele may well be a protective factor 
against AD of greater importance than APOE2, although 
this awaits proof of Thl/E47cs-T allele over-representation 
in healthy individuals of other populations. 
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IN 1993, Allen Roses and colleagues identified APOE4, the 
allele encoding apolipoprotein E4 (apoE4), as a susceptibility 
gene for late-onset Alzheimer's disease (AD) (14,86). Thus, 
subjects homozygous for APOE4 (APOE4/4) are at much 
greater risk of developing AD at a particular age than subjects 
heterozygous for APOE4 (APOE3/4 or APOE2/4), who in 
turn are at greater risk than those carrying APOE3 or APOE2 
alleles. It also appears that subjects carrying one or two 
APOE2 alleles are "protected" from the disease, such that 
the likelihood of developing AD increases much later in life 
than their APOE3 or APOE4 counterparts (13) (see Fig. 1). 
These findings have now been confirmed by many other work- 
ers [see (7,36,60,81,90,98)] and have triggered a wealth of 
interest in investigating the possible role of apoE in the brain, 
and in particular, its role in the development of AD. 

Alzheimer's disease is characterised phenotypically by a 
progressive cognitive decline, in particular a loss of memory, 
learning, and attention (31,57,68). Pathologically, AD is con- 
firmed postmortem by the presence of neuritic amyloid 
plaques and neurofibrilliary tangles (11,84,92). Histological 
and brain imaging studies reveal significant brain shrinkage, 
particularly of temporal lobe structures (43,69,91), and this 
corresponds with morphological changes also reported in this 
region, including decreases in dendritic arborisation (17) and 
loss of neurones and synapses (27,54,109). Regions that show 



Pharmacology Biochemistry and Behavior, Vol. 56, No. 4, pp. 675-685, 1997 

Copyright © 1997 Elsevier Science Inc. 
Printed in the USA. All rights reserved 
0091-3057/97 $17.00 + .00 



the most significant pathology include the temporal lobes and 
the frontal, parietal, and posterior cingulate cortices (9,15), 
areas associated with cognitive function. It now appears that 
loss of synapses and neurofibrillary tangles correlate best with 
the progressive cognitive decline (3,9,15,85,109), indicating 
that neuronal dysfunction is likely to be one of the earliest 
pathologic mediators of the symptomatology. Thus, research 
on apoE has focussed largely on investigating its effects on 
neuronal function and dysfunction, using both in vitro systems 
and in vivo models. In this review, we will discuss the knowl- 
edge to date and consider the evidence suggesting apoE as 
both a physiological mediator of neuronal function and a key 
mediator of pathologic development of AD. 

CNS DISTRIBUTION 

Many of the early investigations into the role of apoE in 
the brain and AD were considered in the context of the well- 
established functions of apoE within the periphery [see 
(47,114)], and it does appear that there are some consistent 
parallels, although interesting differences are emerging. In 
particular, the cellular distribution of apoE within the brain is 
opening new ideas as to the possible function of this ubiquitous 
protein. Thus, plasma apoE originates from many peripheral 
organs, particularly the liver, and is known classically to be 
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AGE AT ONSET 

FIG. 1. Graph to illustrate the influence of the APOE genotype on 
age of onset for Alzheimer's disease. Onset curves were determined 
by Kaplan-Mayer product limit distributions. An indication of the 
influence of the APOE genotype on AD onset may be seen by compar- 
ing incidence rates at age 75; thus, approximately 80% of subjects 
homozygous for APOE4 (4/4) were diagnosed with AD, compared 
with about 40% of subjects homozygous for APOE3 (3/3). The inci- 
dence of AD in patients carrying an APOE2 allele is extremely rare, 
particularly in 2/2 and 2/3 subjects. [Figure reproduced from reference 
(82) with permission of the author.] 



involved in the transport of lipoproteins containing cholesterol 
and lipids to various cells. This function is mediated by the 
binding of apoE to two different receptors: the low density 
lipoprotein (LDL) receptor and the LDL receptor-related 
protein/a2-macroglobulin receptor (LRP). The brain is the 
second most important site for apoE biosynthesis, and for 
many years apoE has been considered as an astrocytic protein 
(47,95), with a key role in the transport of cholesterol. Recent 
studies now reveal that, in addition to its expression in 
astrocytes, apoE can be expressed in neurones, although there 
appear to be considerable species differences. Thus, in rats and 
mice, apoE-like immunoreactivity does not normally appear to 
be expressed in neurones, but rather in microglia and astro- 
cytes. However, in primates (88) and humans (28,29), apoE- 
like immunoreactivity is found in both astrocytes and neu- 
rones. In humans, apoE localises to neurones both in normal 
brain and in AD brain. In AD brain, apoE can be detected 
in neurones both with and without neurofibrillary tangles (29), 
and it has recently been shown that the intensity of apoE- 
like immunoreactivity appears to occur in areas of the brain 
showing the most severe pathology (Einstein et al., unpubl. 
data). At the subcellular level, apoE has been demonstrated 
in several compartments, including the cytoplasm (28,29). 

The difference between species in the cellular localisation 
of apoE is also confirmed by studies in transgenic mice, which 
express the different human apoE isoforms in the absence of 
the native apoE [(117); see Transgenic Research section be- 
low]. In mice expressing either apoE2, apoE3, or apoE4, apoE 
can be detected in both neurones and glia, a pattern of immu- 
noreactivity similar to that for the primate. It is still unclear 
whether apoE is made in neurones rather than being simply 



transported into neurones from another synthetic site. Indeed, 
brain neurones do express LRP receptors (10,79,110,115) as 
well as another member of the LDL receptor family, the very 
low density lipoprotein receptor (VLDLR) (41), and it is sug- 
gested that these receptors may mediate the uptake of apoE 
into neurones. The role of the LDL receptor in mediating 
neuronal uptake of apoE is unclear, as the neuronal expression 
of the LDL receptor remains controversial (79,108). To date, 
there has been no evidence from in situ hybridisation tech- 
niques for the localisation of mRNA for apoE in neurones 
(16). 

CELLULAR AND MOLECULAR INTERACTIONS 

Much of the difficulty in assessing the role of apoE in 
AD stems from the diversity of interactions that the protein 
demonstrates within the central nervous system (CNS). In the 
periphery, the story is fairly clear: apoE mediates the transport 
of lipid between and within cells [see (47)]. However, even 
this role has yet to be definitively proven in the CNS. There 
are practicalreasons why it has been difficult to determine 
the function of apoE in the CNS: the protein is only weakly 
expressed in rodent brain under normal physiological condi- 
tions, so its distribution and site of synthesis, as described in 
the preceding section, remain controversial. Secondly, in vitro 
studies are hampered by the inherent stickiness of the protein, 
which may favour interactions that do not occur under physio- 
logical conditions in vivo. 

Interaction with Microtubule Associated Proteins 

ApoE3 is able to form sodium dodecyl sulfate-resistant 
complexes with two microtubule associated proteins, tau and 
MAP2c, in vitro; apoE4 is unable to form these stable com- 
plexes (106). MAP2c and tau promote the binding and assem- 
bly of microtubule structures, and thus are intrinsic to the 
stabilisation of the neuronal cytoskeleton (23). The formation 
of paired helical filaments and neurofibrillary tangles is a hall- 
mark of AD; therefore, it has been proposed that the increased 
risk of AD in people carrying the APOE4 allele may arise 
from an inability to maintain microtubule integrity (104). 

A feature of the cytoskeleton in the AD brain is the pres- 
ence of hyperphosphorylated tau. It was originally proposed 
that hyperphosphorylation would occur because tau was not 
protected by binding of apoE, and that less protection would 
be afforded by apoE4 than by apoE3, as discussed above. 
However, formation of the apoE-tau complex is also depen- 
dent upon the phosphorylation state of tau, because phosphor- 
ylation of a site within the microtubule binding domain of tau 
has been shown to prevent binding of apoE (37). Thus, one 
might envisage that loss of apoE binding and phosphorylation 
of tau would conspire to cause a rapid breakdown of the 
neuronal cytoskeleton. 

In addition to the formation of a stable complex between 
apoE and tau, a reversible complex, which may be formed 
equally by apoE3 and apoE4, has been shown (103). If the 
formation of this complex is more important in vivo, then a 
rationale for the increased risk of AD conveyed by the APOE4 
gene is diminished; however, the difference between apoE4 
and apoE3 could arise from a difference in subcellular localisa- 
tion rather than any isoform-specific interaction. 

Both the stable and reversible interactions of apoE with the 
cytoskeleton have been shown in vitro, so it will be necessary to 
demonstrate that either one or both occur under physiological 
conditions. In rodents, there is an apparent disparity between 
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the level of apoE expression in neurones compared with that 
for tau or MAP2c, which would argue against a critical role 
for apoE in cytoskeletal maintenance. However, results from 
the in vitro studies show that the formation of apoE-tau com- 
plexes can be detected at very low concentrations of apoE 
[10~ 9 M (106)], and, furthermore, evidence seems to point to 
higher levels of apoE in human neurones (see the next sec- 
tion). The relevance of apoE-tau interactions to normal physi- 
ology, or indeed AD, is far from clear at present. In addition 
to experiments to study apoE function in vivo, it will be impor- 
tant to confirm the correlation between buildup of neurofibril- 
lary tangles and APOE4 genotype, as has been described by 
Ohm et al. (65). 

Regulation of Neurite Outgrowth 

The first study to demonstrate a difference in the activity 
of different apoE isoforms on intact neurones was reported 
by Nathan et al. (61). Treatment of rabbit dorsal root ganglion 
neurones with apoE3 resulted in an increase in neurite exten- 
sion and a decrease in neurite branching. In contrast, apoE4 
treatment reduced both neurite extension and neurite 
branching. Furthermore, the effects of apoE4 predominated 
over apoE3 when combinations of the two isoforms were 
employed (61). These effects of the different apoE isoforms 
were dependent upon the presence of very low density lipopro- 
tein (VLDL), since neither apoE3 nor apoE4 had any effect 
alone, indicating that apoE needs to be presented to neurones 
as a lipoprotein complex to elicit these effects. These effects 
were subsequently reproduced with mouse neuroblastoma 
(Neuro-2a) cells stably expressing human apoE3 or apoE4 
(5). The differential effects of the apoE isoforms were again 
dependent on the presence of VLDL and were blocked by 
agents preventing the uptake of lipoproteins by LRP (5,34). 

Further studies of Neuro-2a cells treated with the different 
apoE isoforms demonstrated that apoE4 treatment resulted 
in production of fewer microtubules and a greatly reduced 
ratio of polymerised to monomeric tubulin than treatment 
with apoE3 (62). It remains to be established whether apoE 
exerts its effects on the neuronal microtubular cytoskeleton 
by a direct or an indirect mechanism. However, the demonstra- 
tion of a differential direct interaction between apoE and the 
microtubule associated proteins MAP2 and tau (106) does 
define a novel mechanistic basis for the differences observed 
in cell cultures treated with the different isoforms. 

An alternative mechanism by which apoE might affect 
neuronal growth and plasticity is interaction with extracellular 
matrix proteins. It has been shown that apoE can interact 
with extracellular matrix proteins such as laminin (38), and 
experiments in which neurones were cultured on substrates 
of either laminin alone or laminin plus apoE showed increased 
cell number and neurite outgrowth in the presence of ApoE. 
However, unlike the isoform-specific effects of apoE in the 
Nathan (61) study, the effect in the presence of laminin was 
identical whether apoE3 or apoE4 was used. 

With respect to the association of apoE4 with late onset 
Alzheimer's disease, the apoE4 isoform may contribute to 
development of the disease by inhibiting neuronal remod- 
elling, as reflected by an inhibition of neurite extension in 
culture. In contrast, the enhancement of neurite extension in 
culture by apoE3 may reflect the ability of this isoform to 
promote remodelling of neuronal processes to compensate for 
the dendritic pruning observed in ageing and Alzheimer's 
disease. 



Lipid Transport 

The accepted role for apoE in the periphery is the media- 
tion of lipid transport (47). ApoE3 and apoE2 tend to associate 
with high density lipoprotein (HDL) particles, whereas apoE4 
tends to bind VLDL particles. The apoE-phospholipid com- 
plexes bind lipoprotein receptors and are internalised. ApoE3 
and apoE4 bind equally well to the LDL and LRP receptors, 
whereas apoE2 binds with 100-fold lower affinity (reviewed 
by Mahley (47)]. Thus, APOE2 homozygotes show elevated 
levels of plasma cholesterol and have a high risk of atheroscel- 
erosis [see (93) for review]. 

Because interaction of apoE with the LDL and LRP recep- 
tors is identical for the APOE3 and APOE4 isoforms, such 
interaction cannot account for the association of the latter 
with AD. However, the intracellular fate of apoE-Iipid com- 
plexes may differ between apoE3 and apoE4, as discussed 
earlier. Although relatively few studies have examined the 
role of apoE in lipid transport in the CNS, the high lipid 
content of the brain (>50%) and the restricted expression of 
members of the apolipoprotein family would likely favour as 
important a lipid-transport role for apoE in the CNS, as in 
the periphery. Furthermore, experimental models of brain 
injury have shown that lipid transport plays a key role in 
subsequent synaptic regeneration and that this process occurs 
concurrently with increased expression of apoE (see later). 
Aspects of lipid trafficking by apoE may also account for some 
of the in vitro interactions described elsewhere in this review. 
For example, in the studies of the effects of apoE on neurite 
outgrowth described above (61), beta-VLDL was required for 
the effect. 

Interaction with ^-amyloid 

The deposition of dense plaques of p-amyloid is a necessary 
pathologic marker of Alzheimer's disease. However, there is 
considerable diversity of opinion regarding the significance of 
plaque formation to the etiology of the disease. Despite this, 
it has been an important objective to demonstrate or refute 
a link between p-amyloid and apoE. Increased deposition of 
P-amyloid has been observed in the cortex of AD patients 
who carried an APOE4 allelle compared with those who did 
not (20,65,87). However, no increase in the rate of progression 
of the disease was reported (87). It is curious that in Down's 
syndrome, a condition that shows many of the pathologic traits 
of AD, no increase in amyloid load was noted in the brains 
of APOE4 carriers, although in this case, the APOE4 genotype 
was associated with earlier death (48). This inconsistency may 
highlight the problems of extrapolating the rate of progression 
of a disease from pathologic markers identified postmortem. 
Increased deposition of p-amyloid in APOE4 carriers has also 
been documented in patients who have died of head injuries 
or stroke (64). This implies that apoE may have a more general 
role in the response to brain injury rather than a specific one 
relating to AD [see (89)]. However, this would be consistent 
with apoE as a risk factor for AD rather than a cause. 

ApoE has been shown to be associated with p-amyloid 
plaques in the AD brain (Fig. 2), although it should be noted 
that other members of the apolipoprotein family, apoJ and 
apoAl, were also present (30), which is suggestive of a simple 
aggregation of lipophilic proteins rather than a specific interac- 
tion. However, a number of groups have studied the interac- 
tion between apoE and p-amyloid. Strittmatter et al. (107) 
have shown that human apoE from plasma could bind 
P-amyloid under cell-free conditions. They showed that the 
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FIG. 2. Adjacent sections of postmortem temporal cortex tissue from 
an 88-year-old female AD patient. One section was immunostained 
with an antibody to p-amyloid (10 jxg/ml; cat. no. 138143 1, Boeh- 
ringer), and an adjacent section was stained with an antibody to apoE 
(1:8000; cat. no. AB947, Chemicon). Note the close overlap between 
the distribution of the two proteins within the neuritic plaque. (Figure 
provided courtesy of Lorraine Fellows, Glaxo Wellcome, Stevenage.) 



interaction may be dependent upon the oxidation state of 
apoE, since no interaction was seen in oxygen-free medium; 
furthermore, they found that apoE4 bound far more rapidly 
to p-amyloid than did apoE3. They suggested that oxidation 
of apoE could occur under inflammatory conditions around 
plaques in the AD brain, and that apoE might enhance the 
aggregation of amyloid. A similar interaction between apoE 
and p-amyloid was also shown by Chan et al. (12) using recom- 
binant human apoE derived from Escherichia coli. However, 
in this study no difference was observed between apoE3 and 
apoE4. In contrast to both of these studies, experiments in 
COS cells in which C-terminal fragments of APP were overex- 
pressed with apoE showed that apoE actually reduced the 



formation of amyloid aggregates (66). The authors saw no 
difference between apoE3 and apoE4, and concluded that 
apoE might inhibit amyloid fibril formation in vivo in an iso- 
form-independent manner. A more detailed study by Evans 
et al. (18) found that apoE3, at physiological concentrations, 
could inhibit amyloid aggregation, with apoE4 showing a simi- 
lar but less potent effect. Evans et al. proposed that the linkage 
between the APOE4 allele and AD was due to the reduced 
ability of APOE4 homozygotes to suppress amyloid forma- 
tion. Some resolution of the inconsistency of the studies of 
apoE-amyloid interactions may come from a study by LaDu 
et al. (45), in which they showed that purification of human 
plasma apoE, which would denature and delipidate the pro- 
tein, reduced differences between apoE3 and apoE4. Thus, 
those studies that have used purified preparations (12,66) may 
produce results that do not reflect the likely interactions of 
apoE in vivo. 

The conformation of p-amyloid may also be a critical factor 
in the interaction with apoE and may account for some of 
the discrepancies of the in vitro studies. Transformation of 
P-amyloid from an a-helical soluble structure to an aggregated 
p-pleated sheet is the defining step in amyloid fibril formation. 
ApoE binds more readily to the p-pleated sheet form (24). 
Human plasma apoE was found to increase the aggregation 
of p-amyloid in vitro, as might be expected from the studies 
of Strittmatter et al. (105) and Chan et al. (12). Significantly, 
however, apoE was able to effect a conformational change of 
p-amyloid from the a-helical form to the p-pleated sheet form, 
as measured by circular dichroism, under conditions that mini- 
mised aggregation of the peptide (100). Thus, the authors 
concluded that apoE actively enhanced p-amyloid aggrega- 
tion, although they were unable to distinguish between the 
different isoforms of apoE. 

One interesting hypothesis regarding the interaction be- 
tween apoE and p-amyloid is that the former may actively 
scavenge amyloid from the extracellular space and thus pre- 
vent its aggregation. Significantly, active microglia, some of 
which were immunoreactive for apoE, have been observed 
in clusters close to the core of amyloid plaques (112). This 
hypothesis is consistent with recent studies of apoE knockout 
transgenic mice, in which Robertson et al. found early deposi- 
tion of p-amyloid in the brains of the mice, suggestive of a 
deficit in their ability to clear extracellular amyloid (unpubl. 
report: Roses, pers. comm.). On the basis of this result, coupled 
with the increased amyloid load observed in APOE4-carrying 
AD patients, it is possible that possession of the APOE4 allele 
is equivalent to a loss of apoE function; therefore, reduced 
interaction with p-amyloid, such as that observed by Evans 
et al. (18) or Ohman et al. (66), might be more relevant to 
Alzheimer's disease. Further studies in transgenic animal mod- 
els of AD to investigate the effects of different human apoE 
isoforms on amyloid load will be required to resolve these 
issues. 

Contribution to Oxidative Stress 

It has been suggested that the generation of free radicals 
may contribute to the neuronal damage seen in Alzheimer's 
disease (6,67). In a study of apoE-deficient mice, Matthews and 
Beal (56) observed increased levels of nitrotyrosine residues 
in various brain regions, which is indicative of heightened 
oxidative stress mediated by peroxynitrate. They proposed 
that a role for apoE in the CNS might be to scavenge reactive 
oxygen species. ApoE2 contains two cysteine residues at posi- 
tions 112 and 158; the Cys 158 residue is mutated to arginine in 
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apoE3. ApoE4 contains no cysteine residues in these positions, 
being identical to E3 except in mutation of Cys to Arg at the 
112 position. The cysteine residues may be readily oxidised 
and so may mop up free radicals. Thus, apoE2 would show 
greater scavenging ability than apoE3, whilst apoE4 would 
not show this activity and so would be more equivalent to the 
knockout. However, the hypothesis is at odds with theories 
to account for the difference in p-amyloid binding of apoE3 
versus apoE4. In the latter case, the rate of binding has been 
shown to be dependent upon oxidation of apoE, and it is 
proposed that apoE4 is more rapidly oxidised than apoE3 
(105). Whatever the relative abilities of apoE3 and apoE4 to 
neutralise free radicals, the 1:1 stoichiometry of this interac- 
tion and the low levels of apoE, at least in rodent neurones, 
would argue against such a role, although neuronal levels of 
apoE may increase dramatically in tissues undergoing trauma 
(see the next section). If the theory is correct, then one might 
expect to see apoE localised in regions of the cell where free 
radicals are likely to be generated, e.g., peroxisomes and mito- 
chondria. However, more specialised mechanisms, such as the 
enzyme superoxide dismutase, exist to deal with free radicals 
in neurones. This is likely to be an area of considerable future 
interest, and the availability of human apoE isoform transgenic 
animals models is likely to be of benefit in this research. 

IN VIVO STUDIES OF ApoE 

In vivo studies on apoE have been somewhat limited to 
date and are focussed in two areas. The first involves the 
production and study of mice in which the APOE gene has 
been either eliminated or replaced by hAPOE isoforms [the 
mouse normally expresses a single apoE isoform most similar 
in structure to apoE4; see (8)]. The second involves the study 
of apoE expression in models of CNS trauma, notably cerebral 
ischaemia, and CNS remodelling after entorhinal cortex 
lesions. 

Transgenic Research 

A valuable means with which to investigate the function 
of genes and their protein products is the technique of gene 
targetting and the generation of transgenic animals overex- 
pressing, or deficient in, the gene of interest. In 1992, two 
groups, Maeda and coworkers (University of North Carolina, 
Chapel Hill, NC) and Breslow and coworkers (Rockefeller 
University, New York), reported the generation of APOE 
knockout mice by homologous recombination to disrupt the 
APOE gene in an ES cell line (129 or E14 derived) (70,71). 
The resulting chimeras were backcrossed to either 129, B ALB/ 
cJ, or C57/B6 lines. As might be expected, the major pheno- 
type of these animals was high plasma triglyceride and choles- 
terol levels, the latter being approximately five times the level 
of wild-type controls (118). This underlies the development 
of fatty streaks and foam cell deposits in the aortic sinus, 
ultimately leading to complex atherosclerotic lesions virtually 
indistinguishable from those seen in humans (71,94,118). 
These may begin to develop in the knockout mice as early as 
3 months of age (71,118). 

In terms of the effect of APOE gene deletion on CNS 
morphology and function, the most comprehensive analysis 
to date has been reported by Masliah and coworkers (53), 
primarily using the Rockefeller-derived stock (71). They re- 
ported significant reductions in dendritic and synaptic density 
in hippocampal and cortical areas, respectively, measured us- 
ing MAP2 and synaptophysin immunoreactivity. Further 



changes in other cytoskeletal markers, e.g., a- and p-tubulin, 
were also noted (53) and may be seen to be consistent with 
the belief that apoE may be involved in the stabilisation of 
microtubule assembly [see (103,104,106)]. These changes be- 
came apparent as early as 4-8 months of age although, perhaps 
surprisingly, they did not increase in magnitude beyond this 
period (53). The synaptophysin results have not been repli- 
cated in a very preliminary study reported by Gandy et al. 
(19) using mice derived from the same transgenic line as that 
used by Masliah, hence further studies are necessary. More 
recently, Michaelson and colleagues have published reports 
suggesting cholinergic deficits and increased AT8 and ALZ50 
immunoreactivity of apoE knockout brain suggestive of tau 
hyperphosphorylation [(21,25); although see (58)]. These lat- 
ter findings are of interest given the proposal that, in humans, 
to account for the association between the E4 genotype and 
AD, there may be isoform differences in tau binding [see 
(106)]. 

Markers of oxidative stress include increases in lipid peroxi- 
dation and nitration of proteins [e.g., (6)]. Matthews and Beal 
(56), reported that apoE knockout mice show increased pro- 
tein nitration, assessed by measuring 3-nitrotyrosine levels by 
HPLC, indicating formation of peroxynitrate, a highly reactive 
and neurotoxic free radical (4,113). These CNS changes are 
largely non-region-specific and do not extend to markers of 
lipid peroxidation, measured by DHB A assay. It is of obvious 
interest to explore these observations further and to examine 
a potential relationship between protein nitration and any sy- 
naptic disruption, as reported by Masliah and coworkers (53). 

Given these reported ultrastructural and biochemical dif- 
ferences between apoE knockouts and their wild-type con- 
trols, one might predict behavioural/cognitive differences, par- 
ticularly because much of this pathology has been reported 
in limbic regions, e.g., the hippocampus and frontoparietal 
cortex. The only full report, albeit brief, that addresses this 
issue, suggests that knockout mice may show cognitive impair- 
ments, as assessed using a modified water maze paradigm 
designed to measure working memory by two-trial acquisition 
to a hidden island, the location of which is changed daily (25). 
Curiously, despite showing poor learning within days (working 
memory), the knockouts showed good learning across days 
(reference memory), the magnitude of which seemed superior 
to controls. We are presently addressing a number of these 
issues, including the behavioural/cognitive differences be- 
tween apoE knockouts and their wild-type littermate controls. 
Initial studies have failed to reveal any differences (see Fig. 
3), although it must be noted that these experiments have 
only been conducted in 2-month-old mice. The data reported 
in the studies of Masliah et al. (53), Gordon et al. (25), and 
Genis et al. (21) were in mice at least 4-8 months of age. We 
are presently allowing a separate colony of these animals to 
age before further behavioural/ultrastructural studies are un- 
dertaken. 

An important point to note about transgenic research, par- 
ticularly in behavioural analyses (but presumably in morpho- 
logical/biochemical studies as well), is the nature of the genetic 
background of the transgenic animal under study. Most gene 
targetting is carried out in ES stem cells derived from the 
mouse strain 129. Because of inherent breeding difficulties 
with the 129 strain, frequently these ES cells are introduced 
into a second strain, e.g., C57/B6. The resulting offspring are 
thus a complex mix of background genes from either strain, 
in addition to the gene target. It is now becoming evident that 
mouse strains may show wide variations in motoric/cognitive 
functions, with 129 and C57/B6 strains being extreme examples 
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FIG. 3. Effect of APOE gene deletion on mouse cognitive behaviour assessed using the Morris water maze. Mice (n = 12 per group) were 
trained to locate a visible platform of variable location (cued task) followed by a hidden platform of fixed location (place task). This was 
followed by a 60-s probe task, in which the platform was removed and the spatial bias shown by the animal for the island quadrant (filled bar) 
was recorded. Note that there was no difference between APOE knockout (ApoE-, ■) and wild-type littermate controls (ApoE+, •), both 
of which showed similar cued and place learning. The APOE genotype was subsequently confirmed by measurement of plasma cholesterol levels 
(ApoE-*-, 2.6 mmol/litre; ApoE-, 12.6 mmol/Utre). (Figure adapted from data provided by Dr. Richard Anderson, Glaxo Wellcome, Stevenage.) 



(22,116). Consequently, it has recently been suggested that 
inadequate backcrossing may lead to transgenics whose re- 
ported phenotypical differences may be more attributable to 
background genes rather than to the gene target (22). It is 
therefore imperative that, in transgenic research, variability 
of genetic background is minimised by repeated backcrossing, 
use of littermate wild-type controls, or use of rescue experi- 
ments as a control where the phenotypical difference between 
a transgenic and its control is eliminated by replacement of 
the missing functional protein (22). In our present studies with 
the apoE knockout mice, we are using animals backcrossed 
at least six times to the C57/B6 strain, as well as wild-type 
littermate controls. 

A further issue of potential importance that has not been 
addressed to date in the apoE transgenic studies is the develop- 
ment of atherosclerotic plaques in the knockout animals and 
the potential for compromised cardiovascular function that 
may contribute to the CNS pathology. There are a few points 
to note. First, longitudinal studies are unable to dissociate 
these events, as each may develop in mice as early as 3^4 
months of age (53,118). Heterozygous apoE-deficient mice, 
which have half the normal apoE gene dosage and develop less 
extensive atherosclerosic lesions, show little or no evidence of 
CNS pathology (53). Finally, the accelerated atherosclerosis 
results in increased plasma lipid peroxidation and oxidative 
stress (32). It remains to be seen whether this could lead 
to the CNS changes reported by Matthews and Beal (56). 
Therefore, there is presently no evidence to eliminate this 
potential association. The ideal tool for this research would 
be an animal with normalised apoE function in the periphery 
but still deficient in the CNS. An elegant approach to this has 
been reported by Linton et al. (46), who described a technique 
in which apoE knockouts received bone marrow transplanta- 



tion of macrophage cells (which synthesise apoE) from wild- 
type mice. These animals demonstrated virtually complete 
absence of atherosclerosis and normalised blood cholesterol 
levels (46), yet presumably retained no central apoE ex- 
pression. 

Clearly of more fundamental interest to apoE research in 
relation to AD is the generation of mice expressing the differ- 
ent human apoE isoforms in the absence of endogenous mouse 
apoE expression. Such transgenics have recently been pro- 
duced (117). A number of lines of varying apoE gene copy 
number have been made, and in each (except two apoE2 lines) 
the elevated plasma cholesterol seen in the knockout animals 
has been substantially corrected, confirming functional expres- 
sion of the transgene. In these animals, the apoE gene is under 
the control of the human promoter, since constructs used 
to make the transgenics consisted of human genomic DNA 
fragments containing regulatory sequences flanking both ends 
of the apoE gene (117). Perhaps a feature of such human 
regulatory sequences is the neuronal localisation of apoE in 
the CNS of these transgenics [similar to human and primate; 
see (28,29,88)] compared with normal mice, where apoE locali- 
sation is largely restricted to astrocytes. Clearly, these animals 
will represent useful tools with which to study apoE isoform 
differences in basic brain biology and function, and in models 
of CNS injury and repair. This work may also raise important 
species differences between the expression patterns of apoE, 
and suggest that nontransgenic rodents may be inappropriate 
models for the study of the role of apoE in humans. 

ApoE in CNS Injury and Repair 

It is the work of Judes Poirier and colleagues that has made 
the greatest contribution to the study of apoE and CNS repair. 
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FIG. 4. Schematic representation of hypothesised cholesterol-phos- 
pholipid recycling process in the dentate gyrus region following an 
entorhinal cortex lesion (based on the work of Poirier and coworkers; 
see the section on ApoE in CNS Injury and Repair). Degenerating 
nerve terminals are internalised and degraded within astrocytes (1). 
Cholesterol is liberated and may combine with newly synthesised 
apoE to form apoE-cholesterol-lipoprotein complexes (2). These 
apoE complexes are then directed into the circulation, after which 
they may be taken up into brain cells requiring lipids, through uptake 
via the LDL, or LRP (or both) receptor pathway (3). Once internali- 
sed, the cholesterol is released (4) for dendritic growth or synaptogene- 
sis. As a result of this internalisation, neuronal cholesterol synthesis 
is repressed via reduced synthesis of HMG-CoA reductase (7). E, 
apoE; PL, phospholipids; CE, cholesterol esters; FC, free cholesterol. 
[Figure reproduced from reference (72) with permission of the author.] 



This work has used a model system in rats of entorhinal cortex 
lesions (ECL) that produce a well documented synaptic loss 
in the denervated dentate gyrus region followed by compensa- 
tory resprouting from local septohippocampal neurones and 
contralateral entorhinal cortex afferents (51,55,101,102). A 
series of experiments initially revealed a marked elevation 
of several RNA transcripts following ECL, one identified as 
GFAP (77) and another encoding a 35-kDa polypeptide subse- 
quently identified as apoE (75). In situ hybridisation showed 
the site of apoE synthesis to be reactive astrocytes (GFAP 
double labelling) in the denervated dentate gyrus (76). Further 
experiments studied the relationship between apoE expres- 
sion and cholesterol uptake/biosynthesis and eventually led 
to the scheme summarised in Fig. 4. Essentially, cholesterol 
released during terminal degeneration following lesioning is 
taken up by astrocytes thatsynthesise apoE. It is hypothesised 
that stable apoE-cholesterol-lipoprotein complexes are 
formed, released, and targeted to local neurones undergoing 
axonal growth or synaptogenesis by uptake via apoE receptors 
[LRP, LDL, or VLDL; see (10,41,79,108,115)], although Poir- 
ier has shown that the LDL receptor is upregulated in response 
to the lesion (72,73). Thus, it is suggested that apoE serves as 
a transporter, carrying lipids between astrocytes and sprouting 
neuronal targets. We have recently undertaken similar experi- 
ments in nontransgenic C57/B6 mice and found the expression 
of apoE to be limited to hippocampal regions undergoing 
synaptogenesis and to precede this process (33), consistent 
with the work of Poirier and coworkers. 

Considering this fundamentally important role for apoE, 
it is not surprising that processes of synaptic repair appear 
to be compromised in apoE-deficient mice after ECL (52). 



However, following sciatic nerve crush, apoE levels increase 
dramatically in genetically unmodified animals (40,47), but 
recovery still persists in apoE knockouts and indeed is indistin- 
guishable from controls (78). It is believed that in the periph- 
ery, unlike the CNS, alternative apolipoproteins may substi- 
tute for apoE, e.g., apoD, apoA-I, and apoA-IV. 

A number of studies have examined the pattern of apoE 
expression following transient cerebral ischaemia in the rat 
and the gerbil (26,35,42,44). Each reports a significant increase 
in apoE immunoreactivity, notably in neurones within the 
CA1/CA2 field, the cell population particularly vulnerable to 
neurodegeneration. In situ hybridisation shows a change in 
apoE mRNA to parallel that of GFAP, suggesting that apoE 
is synthesised by astrocytes (2). It remains contentious as to 
whether neurones actually synthesise apoE or, alternatively, 
if apoE is exclusively derived from astrocytes and taken up 
into neurones via LDL or LRP receptors (see earlier sections). 
Interestingly, the peak of apoE expression coincided with the 
period of cell death, leading to the question of whether apoE 
is expressed in dying CA1/2 neurones as a response that en- 
hances survival, or whether it contributes to the cell death. 

If it is to be believed that apoE has a more general role 
in the CNS response to injury, then it might be predicted that 
the APOE genotype may influence the onset or course of 
other neurodegenerative diseases or acute CNS trauma. There 
seems little evidence to suggest that the APOE genotype is 
linked to Parkinson's disease (49), amyotrophic lateral sclero- 
sis [motorneurone disease (59)], Down's syndrome (39,48,50), 
or hippocampal sclerosis dementia (111). However, prelimi- 
nary reports have suggested an association between the 
APOE4 genotype and poorer outcome following intracerebral 
haemorrhage, as assessed using the Barthel diagnostic index 
(1). Additionally, the APOE4 allele may also serve as a risk 
factor in recovery from head trauma (64,99). Furthermore, 
recovery following coronary bypass surgery is influenced by 
a variety of factors, including APOE genotype (63). This is 
clearly an area of considerable interest, and it is likely that in 
the future a number of similar epidemiological studies should 
help to clarify the issue. 

ApoE and the Cholinergic System 

A possible association between ApoE genotype and cholin- 
ergic dysfunction in AD has been proposed, based on postmor- 
tem examination of brain tissue. Thus, with increasing APOE4 
allele copy number, there is a reduction in hippocampal and 
cortical cholinacetyltransferase (ChAT) levels and reduced 
numbers of acetylcholinesterase-positive neurones in choliner- 
gic cell body regions, including the diagonal band of Broca and 
nucleus basalis of Meynert (72,74). Riekkinen and coworkers 
have also essentially replicated these findings (96). It has been 
proposed that this apparent sensitivity of the cholinergic sys- 
tem to the APOE genotype may be linked to its reliance on 
phospholipid metabolism (72,74). It is interesting to note that 
Gordon et al. (25) reported that mice deficient in apoE to 
show evidence for a central cholinergic hypofunction, although 
we have failed to replicate this observation (Anderson et al., 
unpubl. work). 

The cholinesterase inhibitor tacrine (Cognex®) is presently 
the only drug available for the treatment of AD in the USA. 
There have been a number of clinical studies to evaluate its 
efficacy in AD, and each suggests considerable variation in 
patient response to this drug. Recently, Poirier and colleagues 
have suggested that part of this variability may be due to the 
APOE genotype, with patients carrying one or two copies of 
the E4 allele being poor responders (74), the explanation being 
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that these individuals show greater atrophy of the cholinergic 
system. An alternative explanation is misdiagnosis: on post- 
mortem examination, approximately 30% of non-E4 patients 
diagnosed with probable AD show incorrect diagnosis, for 
they display pathology typical of other dementias such as 
Pick's disease, hippocampal sclerosis dementia (83). In con- 
trast, E4 carriers diagnosed with probable AD show about 
98-100% accuracy of diagnosis (83). Hence, tacrine may be 
more beneficial in the non-AD cases. Future studies with 
newer cholinomimetics will undoubtedly establish the ro- 
bustness of these associations. In any case, these data highlight 
an issue that is likely to become increasingly important in 
clinical trial design: that is, an individual's genotype may in- 
fluence treatment response and in future may determine the 
treatment strategy being used. 

CONCLUSION 

The fundamental question concerning the role of apoE in 
the CNS is whether it serves a physiological and/or pathologic 
role. Studies of the expression of apoE in disease tissue, includ- 
ing AD brain, necessarily address a pathologic role for the 
protein; similarly, studies of apoE in CNS injury and repair 
provide evidence that the protein is recruited to stressed or 
dying neurones, either to support regeneration or to oversee 
their demise. In this respect, epidemiological studies designed 
to clarify the influence of the APOE genotype on other neuro- 
degenerative conditions will be of interest. 

The influence of the APOE genotype on Alzheimer's dis- 
ease is clearly demonstrated. However, the question remains, 
does apoE contribute to the onset and early stages of disease, 
or does it simply affect the progression of the disease? Associa- 



tion of apoE with diverse pathologic states would seem to 
point to the latter. However, recent data from PET studies 
showing abnormal glucose metabolism in the cortex of APOE4 
homozygotes up to 20 years before they develop symptoms 
of disease (80) are the first indication that the protein may be 
involved early in the disease process. This is supported by 
observations that in a nondemented aged population, persons 
carrying at least one APOE4 allele had poorer cognitive test 
scores and smaller hippocampal volume (97). 

At the preclinical research level, there is a clear need for 
further study of potential roles for apoE in the CNS, such as 
lipid transport, as an antioxidant, or regulation of neuronal 
cytoskeleton. These lines of investigation are presently quite 
controversial, but will benefit from a better understanding of 
the expression and localisation of apoE in neurones and glia. 
Transgenic studies are still in their infancy, but the availability 
of different lines harbouring human apoE isoforms (117) will 
doubtless prove a valuable means to investigate isoform differ- 
ences in both physiological and pathologic functions of apoE. 

These are early days for apoE research, although there has 
been rapid progress driven largely by the recent estimates 
that approximately 60-70% of all AD cases are linked to the 
APOE4 genotype. There is considerable optimism that an 
understanding of how apoE-isoform differences can modify 
the risk of AD will eventually lead to an effective treatment 
for this devastating disease. 
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Early protective effect of CCR-5 delta 32 heterozygosity on HIV-1 
disease 

progression: relationship with viral load. The SEROCO Study Group. 

Meyer L, Magierowska M, Hubert JB, Rouzioux C, Deveau C, Sanson F, 
Debre P, 

Delfraissy JF, Theodorou I. 

OBJECTIVE: To determine the influence of heterozygosity for the delta 
32 mutant 

CCR-5 allele on HIV-1 disease progression. DESIGN: HIV-1 disease 
progression and 

serum viral load were analysed according to the C-C chemokine receptor 
(CCR)-5 

genotype in 412 Caucasian patients (319 men and 93 women) with a known 
date of 

seroconversion, who were enrolled in the SEROCO cohort (median follow- 
up, 74 

months). RESULTS : The frequency of heterozygosity for the mutant allele 
was 17% 

and did not differ according to sex or risk factor of HIV infection. 
Heterozygotes were significantly less likely than patients with two 
functional 

alleles to have symptomatic primary infection. Their serum viral load 
was lower 

during the 6- to 24 -month plateau phase after seroconversion. This 
difference 

persisted afterwards, although the rate of decline in CD4+ cells was 
similar. 

Kaplan-Meier survival curves showed slower progression to clinical AIDS 
in 

heterozygotes during the first 7 years following infection (P < 0.02), 
the two 

curves tending to join thereafter (overall log-rank test, P = 0.17). 
However, 

the interaction term with time did not reach significance in a Cox 
model . The 

overall relative risk of progression was 0.67 (95% confidence interval, 
0.38-1.18) and was not influenced by adjustment for age at 
seroconversion or 

symptomatic primary infection. After adjustment for early viral load 
the 

relative risk was 0.83. Pneumocystis carinii pneumonia and 
toxoplasmosis were 

less likely to be the first AIDS-def ining illness in heterozygotes than 
in the 

other patients (0 versus 24.7% of AIDS cases, P = 0.04), despite 
similar 

management. CONCLUSION: Deletion of one CCR-5 gene allele appears to 
protect 

against HIV-1 disease progression, mainly during the early years of the 
infection. Heterozygosity for the deletion leads to persistently lower 
viral 

load, and also seems to protect against some opportunistic infections. 
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Evidence for a Major Quantitative Trait Locus on 
Chromosome 17q21 Affecting Low-Density Lipoprotein Peak 

Particle Diameter 

Yohan Bosse, MSc; Louis Perusse, PhD; Jean-Pierre Despres, PhD; Benoit Lamarche, PhD; 
Yvon C. Chagnon, PhD; Treva Rice, PhD; D.C. Rao, PhD; 
Claude Bouchard, PhD; Marie-Claude Vohl, PhD 

Background — Several lines of evidence suggest that small dense LDL particles are associated with the risk of coronary 
heart disease. Heritability and segregation studies suggest that LDL particle size is characterized by a large genetic 
contribution and the presence of a putative major genetic locus. However, association and linkage analyses have thus 
far been inconclusive in identifying the underlying gene(s). 

Methods and Results — An autosomal genome-wide scan for LDL peak particle diameter (LDL-PPD) was performed in the 
Quebec Family Study. A total of 442 markers were genotyped, with an average intermarker distance of 7.2 cM. 
LDL-PPD was measured by gradient gel electrophoresis in 681 subjects from 236 nuclear families. Linkage was tested 
by both sib-pair- based and variance components-based linkage methods. The strongest evidence of linkage was found 
on chromosome 17q21 .33 at marker D17S1301, with an LOD score of 6.76 by the variance-components method for the 
phenotype adjusted for age, body mass index, and triglyceride levels. Similar results were obtained with the sib-pair 
method (P<0.0001). Other chromosomal regions harboring markers with highly suggestive evidence of linkage 
(P<0.0023; LOD > 1.75) include lp31, 2q33.2, 4pl5.2, 5ql2.3,and 14q31. Several candidate genes are localized under 
the peak linkages, including apolipoprotein H on chromosome 17q, the apolipoprotein E receptor 2, and members of the 
phospholipase A 2 family on chromosome lp as well as HMG-CoA reductase on chromosome 5q. 

Conclusions — This genome- wide scan for LDL-PPD indicates the presence of a major quantitative trait locus located on 
chromosome 17q and others interesting loci influencing the phenotype. (Circulation. 2003;107:2161-2368.) 

Key Words: genome ■ lipoproteins ■ genetics ■ genes 



A number of case-control as well as prospective studies 
reveal an increased risk of coronary heart disease in 
patients with small dense LDL compared with those having 
larger, more buoyant LDL particles. 1 Heritability studies 
based on twins suggested that approximately one third to one 
half of the variation in the LDL peak particle size can be 
attributed to genetic influences. 2 ' 3 Complex segregation anal- 
yses of small dense LDL phenotypes have been performed 
with data from different types of family structures, different 
criteria for proband ascertainment, and the use of different 
techniques to characterize LDL heterogeneity. 4 - 9 Neverthe- 
less, the model providing the best fit to the data included 
either a dominant, a recessive, or an undetermined mode of 
inheritance for the trait. Furthermore, the allele frequency 
determining the small dense LDL phenotype ranges from 
19% to 42%, with reduced penetrances in young men and 



premenopausal women. However, these studies unanimously 
provided evidence in favor of a gene with a major effect on 
LDL particle phenotypes. 

Association studies with candidate genes have been incon- 
sistent in finding genes associated with small dense LDL, The 
-250G— >A polymorphism within the hepatic lipase pro- 
moter was associated with buoyant LDL particles. 10 How- 
ever, the -5140— »T polymorphism, which is in complete 
linkage disequilibrium with the -250G— >A polymorphism, 11 
showed no effect on LDL particle size. 1213 The apolipopro- 
tein (apo) E genotype was also associated with the small 
dense LDL phenotype. However, some studies have reported 
smaller particles for subjects carrying the E4 allele, 14 - 16 and 
other studies did for subjects carrying the E2 allele. 17 - 1 8 In 
contrast, others have shown that LDL particle size did not 
differ among the apoE genotypes. 19 Additional candidate 
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TABLE 1. Descriptive Statistics of LDL Peak Particle Diameter and Covariates in 
Each of the Sex and Generation Groups 

Fathers Mothers Sons Daughters 

Variables (n=137) (n=194) (n=148) - (n=202) 



LDL-PPD, A 


261.2 


±5.4 


264.3 


±5.2 


262.6± 


4.5 


264.4± 


4.5 


Age, y 


55.5 


±9.2 


55.5 


±12.7 


26.6 ± 


10.0 


28.2± 


10.7 


BMI, kg/m 2 


28.9 


±6.3 


28,7 


±8.0 


26.5± 


6.9 


27.6 ± 


8.7 


Triglyceride, mmol/L 


1.98 


±1.23 


1.79 


±2.40 


1.35± 


0.71 


1.27± 


0.59 



Values are mean±SD. 



genes, including cholesteryl ester transfer protein, 20 microso- 
mal triglyceride transfer protein, 21 cholesterol 7a- 
hydroxylase, 22 apoB-100, 2 - 1 apoC-III, 24 and ACE 16 were in- 
vestigated for potential effects on small dense LDL 
phenotypes. These studies revealed either the absence of an 
association or the presence of an association only in particular 
subgroups. 

Results from linkage studies are equivocal. After linkage of 
small dense LDL with the apoB (the protein moiety of LDL) 
gene locus on chromosome 2 was excluded, 25 ' 26 suggestive 
linkage to the LDL receptor locus on chromosome 19 has 
been reported. 27 * 28 However, subsequent sequencing of the 
entire coding regions of the LDL-receptor gene did not reveal 
any sequence variants, thus weakening the hypothesis that a 
mutant LDL-receptor allele is responsible for the dense LDL 
phenotype. 29 Other candidate loci, including hepatic lipase, 12 
lipoprotein lipase, 30 cholesteryl ester transfer protein, 28 - 31 - 32 
apoAl-CIH-AIV complex, 28 * 32 and the manganese superoxide 
dismutase, 28 - 32 have been shown to be linked with the small 
dense LDL phenotype. Unfortunately, most of these linkages 
have not been replicated. 33 * 34 On the basis of these results, 
Austin et al 34 emphasized the necessity of finding new 
genetic loci, other than those harboring known candidate 
genes, to identify the genes potentially involved in determin- 
ing the small dense LDL phenotype. Genome-wide scans are 
particularly suited for this purpose. Previous genome-wide 
scans have focused on variation in cholesterol concentrations 
of LDL size fractions. Rainwater et al 35 found 2 quantitative 
trait loci (QTLs) on chromosome 3 and 4 with logarithm of 
the odds (LOD) scores >3 for LDL size fraction 3 (LDL-3), 
a fraction that contains small LDL particles. This study 
demonstrates the existence of QTLs affecting the concentra- 
tion of cholesterol within a particular subpopulation of LDL 
but do not provide evidence of QTLs responsible for the size 
of the LDL particle by itself. To the best of our knowledge, 
the only whole-genome scans on LDL particle size were 
performed on 240 individuals ascertained through 18 unre- 
lated familial combined hyperlipidemic probands, 12 Results 
suggest a locus over the hepatic lipase gene on chromosome 
15 with a LOD score of 2.2. Here, we report the results of an 
autosomal genomic scan for LDL peak particle diameter 
(LDL-PPD) measured by gradient gel electrophoresis. 

Methods 

Population 

The Quebec Family Study is an ongoing project composed of 
French-Canadian families that has been described previously. 36 In 
the present study, a total of 681 subjects from 236 nuclear families 



had available data on LDL-PPD. Table 1 presents the characteristics 
of subjects in each of the sex and generation groups. The Laval 
University Medical Ethics Committee approved the study, and all 
subjects provided written informed consent. All the procedures 
followed were in accordance with institutional guidelines. 

Phenotype 

LDL-PPD was measured by gradient gel electrophoresis from 
plasma obtained after a 12-hour fast. Details on the technique have 
been provided previously. 37 

Genotypes 

Genomic DNA was prepared by the proteinase K and phenol/ 
chloroform technique. DNA preparation, polymerase chain reaction 
conditions, and genotyping are described in detail elsewhere. 38 
Genotypes for each marker were typed by use of automatic DNA 
sequencers and the computer software SAGA from LICOR. The 
results were stored in a local dBase IV database, GENEMARK, 
which inspects results for mendelian inheritance incompatibilities 
within nuclear families and extended pedigrees. A total of 335 
microsatellite markers (dinucleotide, trinucleotide, and tetranucleo- 
tide repeats) selected from different sources, but primarily from the 
Marshfield panel version 8a, were available for this genome scan. 
The location of markers on the chromosomes in centimorgans (cM) 
were taken from version 9.0 of the Marshfield Institute map 39 and the 
Location Database map. 40 In addition, 107 polymorphisms in 63 
candidate genes were included. The average intermarker distance for 
the whole set of 442 markers was 7.2 cM. The Genome Database 41 
and the OM1M gene map 42 were used to identify candidate genes. 

Statistical Analyses 

LDL-PPD was adjusted for covariates by use of a stepwise multiple 
regression procedure, retaining only terms that were significant at the 
5% level. Regression parameters were estimated within 6 age-by-sex 
(<30, 30 to 50, and ^50 years; male versus female) groups after 
exclusion of outliers (±4 SD), and residuals were computed for all 
subjects. Residual scores were then standardized to a mean of 0 and 
an SD of 1. LDL-PPDs were adjusted for 3 different sets of 
covariates: (1) age up to the cubic polynomial, (2) age and body mass 
index (BMI), and (3) age, BMI, and triglyceride levels. These 
adjustments gave 3 phenotypes, arbitrary called LDL-PPD 1, LDL- 
PPD2, and LDL-PPD3, respectively. The phenotypes were adjusted 
by use of SAS software (version 8.02). 

The search for linkage between the phenotypes and the genetic 
markers was performed by 2 different approaches. First, linkage was 
tested by the new Haseman-Elston regression-based method, which 
models the trait covariance between sib-pairs, instead of the squared 
sib-pair trait difference used in the original method. It regresses the 
mean-corrected trait cross product on the number of alleles shared 
identical by descent (IBD). Single-point and multipoint estimates of 
alleles shared IBD were generated with the GEN1BD software, and 
linkage was tested with the SIBPAL2 software from the SAGE 4.0 
statistical package (SAGE, 200 1). 43 The maximum number of 
sib-pairs was 352. Linkage was also investigated by the variance 
components- based approach implemented in the quantitative trans- 
mission disequilibrium test (QTDT) computer software 44 By this 
approach, the phenotypic covariance among members of a family is 
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Figure 1. Quantitative transmission disequilibrium test linkage results for all autosomal chromosomes with LDL-PPD phenotypes. LOD 
scores are presented on y axis, and genetic distance is presented on x axis in centimorgans. LDL-PPD1, LDL-PPD2, and LDL-PPD3 
indicate LDL-PPD adjusted for (1) age, (2) age and BMI, and (3) age, BMI, and triglyceride, respectively. 



assumed to result from the additive effects of linkage caused by a 
major locus (a), a residual familial component resulting from 
polygenes (g), and a residual nonshared environmental component 
(e) that represents environmental effects unique to each family 
member. Linkage is tested by contrasting the null hypothesis of no 
linkage (ct b =0) with the alternative hypothesis (oa^O) by use of a 
likelihood ratio test, as described previously. 45 The LOD score was 
computed as ^/(2 loge 10). The interpretation of linkage evidence 
was considered suggestive (P<0.0 1 ; LOD > 1 . 1 8), highly suggestive 
(P<0.0023; LOD a 1.75), or evidence of linkage (P<0.0001; LOD 
>3.0). 46 



Results 

An overview of the variance components- based linkage 
results for the 3 LDL-PPD phenotypes is given in Figure 1 . 
Suggestive evidence of linkages (P<0.01 or LOD scores 
>1.18 for at least one of the phenotypes) is summarized in 
Table 2. The strongest evidence of linkage, which was 
confirmed by both linkage methods, was found on chromo- 
some 17q21.33. As shown in Figure 2, the peak linkages were 
found with marker D17S1301 for LDL-PPD 1 (LOD=4.72), 
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TABLE 2. Results From the Genome Scan: Markers Showing Evidence of Linkage With the LDL-PPD Phenotypes According to the 
Linkage Methods Used 
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TABLE 2. Continued 
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LDL-PPD1, LDL-PPD2, and LDL-PPD3 indicate LDL-PPD adjusted for (1) age, (2) age and BMI, and (3) age, BMI, and triglyceride, respectively. Results from single 
point and multipoint are reported in P values, and results from the variance components- based method are reported in L0D scores. Markers with highly suggestive 
evidence of linkage (P<0.0023 or L0D score >1.75) for the three linkage methods for one of the LDL-PPD phenotypes are indicated in boldface type. 

*ft=0.01 or LOD score >1.18. 

1P<0.0023 or LOD score >1.75. 

1P<0.0001 or LOD score >3.00. 

§Chromosome locations for marker with D number or marker with gene name are from the Location Database map, 40 and the other markers are from the Marshfield 
Institute map. 39 



LDL-PPD2 (LOD=4.70), and LDL-PPD3 (LOD=6.76). 
Marker D17S1290, located 1.6 cM from D17S1301, also 
gave fairly good evidence of linkage for the 3 phenotypes. 

Other chromosomes exhibiting some evidence of linkage 
by the variance components-based method are displayed in 
Figure 1. Highly suggestive evidence of linkages was ob- 
served at lp31 (leptin receptor locus), 2q33.2 (marker 
D2S1384), 4pl5.2 (D4S2397), 5ql2.3 (D5S1501), and 
14q31.1 (D14S53). Markers at the leptin receptor locus and 
markers D5S1501 and D14S53 also provided evidence of 
linkage by the sib-pair method (see Table 2). 

Other markers gave highly suggestive evidence of linkage 
(P<0.0023) with at least one of the linkage methods. For 



instance, marker D16S261 provided evidence of single-point 
linkage with the 3 phenotypes. The marker VWFP1 on 
chromosome 22ql 1.21 provided evidence of single-point and 
multipoint linkage for the 3 phenotypes. Conversely, marker 
D4S1627 yielded highly suggestive evidence of linkage for 
LDL-PPD 1 and LDL-PPD2 with the variance component 
method. D5S1457 at 5pl2 shows highly suggestive evidence 
of linkage in multipoint analysis for the 3 phenotypes and in 
single-point for LDL-PPD3. Finally, several markers pro- 
vided highly suggestive evidence of linkage with LDL-PPD3 
only, including D1S198, D2S434, IRS1 (2q36.3), ADRB2 
(5q31), TNFat (6p2L3), D8S1110, D9S1121, D16S410, and 
ACE (17q23). 
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Figure 2. Quantitative transmission disequilibrium test linkage 
results for chromosome 1 7 with LDL-PPD phenotypes. Genetics 
markers used for linkage are indicated under x axis. Approxi- 
mate location of candidate genes in vicinity of major peak are 
displayed on graph. Dashed horizontal line represents a LOD 
score of 3.00. LDL-PPD1, LDL-PPD2, and LDL-PPD3 indicate 
LDL-PPD adjusted for (1) age, (2) age and BMI, and (3) age, 
BMI, and triglyceride, respectively. 



Discussion 

The primary objective of this study was to identify QTLs 
affecting LDL-PPD variation. The results provide evidence 
for a major locus affecting LDL-PPD located on chromosome 
17q21. Interestingly, none of the candidate genes located in 
the area of this QTL were tested previously. The marker 
D17S1301 located on chromosome 17q21 .33 was strongly 
linked with the LDL-PPD, whether adjusted for covariates or 
not. However, the evidence for linkage was stronger when the 
phenotype was adjusted for plasma triglycerides, indicating 
that triglyceride levels may attenuate the penetrance of the 
locus. Marker D17S1290, located 1.6 cM from D17S1301, 
also provided good evidence of linkage (1.34<LOD<2.63). 
The apoH (APOH) gene, also referred to as ^-glycoprotein I, 
is encoded under the peak linkage on 17q21. ApoH is a 
single-chain glycoprotein that exists in plasma both in a free 
form and in combination with lipoprotein particles. It has 
been implicated in several physiological pathways, including 
lipid metabolism, coagulation, and the production of an- 
tiphospholipid antibodies. This apolipoprotein activates li- 
poprotein lipase, 47 and genetic variations in this gene have 
been associated with variation in HDL cholesterol and tri- 
glyceride levels. 48-50 ACE is also located in this genomic 
region. This enzyme cleaves the final intravascular step, 
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TABLE 3. Candidate Genes Within Chromosomic Regions Linked to LDL-PPD 



LOD Score 



Chromosome 



Region 
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LDL-PPD2 


LDL-PPD3 
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2.04 


2.56 


PLA2G2D, PLA2G2A, PLA2G5, AP0ER2, FABP3, CPT2, ABCD3, HMGCS2 


2q33.2 


D2S1384 


1.91 


2.11 


2.27 


PLA2R1 , ABCB11, LRP2, ACADL, FACL3, ABCB6, IRS1, HDLBP 


4p15.2 


D4S2397 


2.14 


2.22 


2.11 


LRPAP1 


5q12.3 


D5S1501 


2.40 


2.34 


2.10 


HMGCR, HMGCS1 


14q31.1 


D14S53 


2.79 


2.65 


1.72 


ABCD4, CYP46 


17q21.32-q21.33 


D17S1290 


1.39 


1.34 


2.63 


ACACA, ABCC3, ACE, APOH, ACOX1 , FASN 




D17S1301 


4.72 


4.70 


6.76 





LDL-PPD1, LDL-PPD2, and LDL-PPD3 indicate LDL-PPD adjusted for (1) age, (2) age and BMI, and (3) age, BMI, and triglyceride, respectively. 



resulting in the vasoactive peptide angiotensin II. Angiotensin 
II has been shown to bind specifically to LDL, 51 which 
produces a modified form of LDL that is taken up by 
macrophages at an enhanced rate, leading to cellular choles- 
terol accumulation. 52 In the present study, the insertion/ 
deletion polymorphism in intron 16 of the ACE gene pro- 
vided evidence of linkage with LDL-PPD1 (LOD=1.46), 
LDL-PPD2 (LOD=L57) and LDL-PPD3 (LOD=2.35). Fig- 
ure 2 shows the approximate location of candidate genes 
surrounding the major peak on chromosome 17. 

Several other chromosomal regions provided highly sug- 
gestive (/ ) <0.0023) evidence of linkage. These regions in- 
clude chromosomes lp31, 5pl2-pl2.3, and 14q31,l, which 
show evidence of linkage with both linkage methods and for 
all LDL-PPD phenotypes. Some promising candidate genes 
are located within these regions. First, the strongest evidence 
of linkage on chromosome lp comes from a marker located 
within the leptin receptor (LEPR) gene. By modulating the 
hypothalamic effects of leptin on food intake and energy 
expenditure, genetic variants in the LEPR may affect energy 
balance and the size of LDL particles as a consequence of 
body fatness alterations. However, adjusting the LDL-PPD 
for BMI did not affect the strength of the linkage. On lp, 3 
members of the phospholipase A 2 (PLA2) gene family are 
present, namely PLA 2 group I ID (PLA2G2D), group V 
(PLA2G5), and group IIA (PLA2G2A). PLA 2 is known to 
hydrolyze the phospholipid monolayers of LDL particles and 
change their physicochemical properties and size. 53 ApoE 
receptor 2 (APOER2) is also located near the locus of 
interest. On chromosome 5, 2 markers (D5S1457 and 
D5S1501) located 20 cM apart provided evidence of linkage 
with LDL-PPD. This region contains the HMG CoA reduc- 
tase (HMGCR), which is the rate-limiting enzyme for cho- 
lesterol synthesis. A list of other potential candidate genes 
within the chromosomal regions linked to the LDL-PPD is 
provided in Table 3. 

Among the panel of markers included in the genome scan, 
few candidate genes for LDL-PPD were present. First, an 
apoB marker gave no evidence of linkage with the phenotype. 
A significant linkage to apoB has been reported in a sib-pair 
linkage analysis of dizygotic female twins, 33 but other linkage 
studies excluded the hypothesis of linkage for the apoB locus 
and LDL size. 25 - 26 * 28 Second, although no linkage was found 
with the lipoprotein lipase locus in the present study or in 2 



others, 28 - 34 a highly significant LOD score of 6.24 was 
obtained in another study of heterozygous lipoprotein lipase- 
deficient families. 30 Third, the apoE gene gave no evidence of 
linkage, as reported previously . 28 - 33 - 34 Finally, consistent with 
3 other studies, 32-34 the LDL receptor also was not linked to 
LDL-PPD in the present study. In contrast, 2 previous results 
linked the LDL receptor locus to LDL subclass in families 
ascertained through probands with the atherogenic lipopro- 
tein phenotype 27 and in families with coronary heart dis- 
ease. 28 However, no amino acid sequence changes in the LDL 
receptor were found in the former study, 27 making it unlikely 
that a mutant allele in the LDL receptor gene was responsible 
for the linkage. 29 In the present study, negative results were 
also obtained with other candidate genes, including paraoxo- 
nase, hormone-sensitive lipase, CD36, and the intestinal fatty 
acid-binding protein. 

In conclusion, the results of this study reveal the presence 
of a major locus located on chromosome 17q2 1 .33 influenc- 
ing LDL-PPD. This finding supports results from a handful of 
segregation analyses indicating the presence of a putative 
major locus for LDL particle size. Evidence of linkage was 
also found on chromosomes lp31, 2q33.2, 4pl5.2, 5ql2.3, 
and 14q31 . 1 . These QTLs harbor a good number of candidate 
genes that have not previously been tested in association 
studies with LDL-PPD. 
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